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ABSTRACT 

We present data and initial results from VLT/X-Shooter emission-line spectroscopy of 96 galaxies selected by long y-ray bursts 
(GRBs) at 0.1 < z < 3.6, the largest sample of GRB host spectra available to date. Most of our GRBs were detected by Swift and 
76% are at 0.5 < z < 2.5 with a median Zmcd ~ 1-6. Based on Balmer and/or forbidden lines of oxygen, nitrogen, and neon, we 
measure systemic redshifts, star formation rates (SFR), visual attenuations (Ay), oxygen abundances (12 + log(0/H)), and emission¬ 
line widths (<t). We study GRB hosts up to z ~ 3.5 and find a strong change in their typical physical properties with redshift. The 
median SFR of our GRB hosts increases from ~ Q.dMeyr^' at z ~ 0.6 up to SFRyaed ~ 15Moyr“* at z ~ 2. A higher 

ratio of [O iii]/[0 ii] at higher redshifts leads to an increasing distance of GRB-selected galaxies to the locus of local galaxies in the 
Baldwin-Phillips-Terlevich diagram. There is weak evidence for a redshift evolution in Ay and cr, with the highest values seen at 
z ~ 1.5 (Ay) or z ~ 2 (cr). Oxygen abundances of the galaxies are distributed between 12 + log(0/H) = 7.9 and 12 -l- log(0/H) = 9.0 
with a median 12 -l- log(0/H)n,ed ~ 8.5. The fraction of GRB-selected galaxies with super-solar metallicities is ~ 20% at z < 1 in 
the adopted metallicity scale. This is significantly less than the fraction of total star formation in similar galaxies, illustrating that 
GRBs are scarce in high metallicity environments. At z ~ 3, sensitivity limits us to probing only the most luminous GRB hosts for 
which we derive metallicities of Z < 0.5 Zq. Together with a high incidence of Z ~ 0.5 Zq galaxies at z ~ 1.5, this indicates that 
a metallicity dependence at low redshift will not be dominant at z ~ 3. Significant correlations exist between the hosts’ physical 
properties. Oxygen abundance, for example, relates to Ay (12 -l- log(0/H) oc 0.17 • Ay), line width (12-1- log(0/H) oc cr°-*), and SFR 
(12 -t- log(0/H) oc S FR^-^). In the last two cases, the normalization of the relations shift to lower metallicities at z > 2 by ~ 0.4 dex. 
These properties of GRB hosts and their evolution with redshift can be understood in a cosmological context of star-forming galaxies 
and a picture in which the hosts’ properties at low redshift are influenced by the tendency of GRBs to avoid the most metal-rich 
environments. 

Key words. Gamma-ray burst: general. Galaxies: evolution. Galaxies: high-redshift. Galaxies: star formation 


1. Introduction 

The extreme luminosities of long y-ray bursts (GRBs) and their 
afterglows (e.g., Gehrels et al. 2009 [ Kumar & Zhang 2015[ 
for recent reviews) over the full electromagnetic spectrum make 
them powerful messengers from the early Universe, a potential 
that was identihed very early on (e.g., Wijers et al.|1998[|Lamb 
|& Reich^|2000| l. Coupled with the association of long GRBs 

* Based on observations at ESO, Program IDs: 084.A-0260, 084.A- 
0303, 085.A-0009, 086.B-0954, 086.A-0533, 086.A-0874, 087.A- 
0055, 087.A-0451, 087.B-0737, 088.A-0051, 088.A-0644, 089.A- 


with core-collapse supernovae of type Ic (e.g., |Galama et al.| 
1998[|Hjorthet al. 2003 [ Plan et al. 2006 1 , and thus the formation 
of some massive stars, GRBs provide means of studying star¬ 
forming galaxies at the highest redshifts (e.g., [Chary et al.|2007[ 
Chen et al.|2009[[Tanvir et al.|2012 l, faintest luminosities (e.g., 
Vreeswijk et al.|2001 Trend et al.||2012| l, as well as the global 
star formation rate (SFR) density up to the epoch of re-ionization 


(e.g., Kistler et al.||2009| Robertson & Ellis||2012| 

Elliott et al. 

2012jl. In particular the information that is obtainec 

through af- 


terglow spectroscopy provides insights of unprecedented detail 
into the chemical composition of high-redshift galaxies (e.g.. 


0067, 089.A-0120, 089.D-0256, 089.A-0868, 090.A-0088, 090.A- 

Fynbo et al. 

2006 

Prochaska et aL|2009| Thone et aL|2013| Cue- 

0760,090.A-0825,091.A-0342,091.A-0703,091.A-0877,091.C-0934, 

chiara et al. 

2015 

1 . 


^ The reduced spectra are only available at the CDS via anony¬ 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u- 
strasbg.fr/viz-bin/qcat? J/A-I-A/5 81/A 125 


Long y-ray bursts, however, are rare objects (e.g., Podsi- 
jadlowsld et al.|[20()4[ [Guetta & Della VallepOOTl ). The rate at 
which GRBs are accurately localized through satellites is around 
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two per week. Number statistics of GRBs are thus orders of 
magnitude lower than those of other prominent probes of the 


early Universe, notably Lyman-break galaxies (LBGs, e.g., Stei 
del et al.|1996|), and Damped-Lyman-a Absorbers (DLAs, e.g., 


Wolfe et al. 112005 1 along the sightline to Quasi-Stellar Objects 


(QSOs). Quite similarly, the GRB explosion represents a very 
rare endpoint of stellar evolution, and could thus be subject to en¬ 
vironmental factors enhancing or quenching the GRB rate with 
respect to star formation. For instance, metallicity often has been 
discussed as being likely to influence GRBs (e.g., |Stanek et al.j 
2006 [ Modjaz et al.|2008 1. Variations in the initial-mass function 


(IMF) or the fraction of massive stars in tight binaries could also 
affect the GRB rate ( [Kelly et al.|2014) l. 

A physical understanding of GRBs and the conditions in 
which they form is thus necessary to put observations of high- 
redshift GRBs into a cosmological context. Arguably, one of the 
key pieces to the puzzle of understanding GRB progenitors and 
their role in probing distant star formation lies in the nature of 


their host galaxies (e.g., Levesque 

20141 for a review). 

Seminal sample studies (e.g.. 

Le Eloc’h et al. 

20031 Chris- 

[tensen et al.[2004|[Tanvir et al. 2004| Levesque et a 

i.|2010a|l cul- 

minating in the compilation of GRB host properties of Savaglio 


-mass. 


star-forming galaxies at z < 1. The host’s UV emission is highly 
concentrated around the GRB position and their morphologies 
appear more irregular than the hosts of core-collapse supernovae 
( jFruchter et al.|200^|Svensson et al.|2010| ). 

In line with this picture were Lya properties ( jFynbo et al.j 
2003 1 Milvang-Jensen et al.||2012l and early measurements of 


sub-solar GRB host metallicities (e.g., jProchaska et al. 2004[ 


[Gorosabel et al.j 2005 [ jSollerman et al.||2005| l. These data sup¬ 

ported a strong dependence of GRB formation with metallicity 
(e.g., [Hirschi et aL||2005[ Yoon et ^|2006|) as would be theo¬ 


retically expected in the collapsar model ( Woosley 1 1993 1 [Mac- 

Fadyen & Woosley|1999 1. 


More recently, however, this previously quite uniform pic¬ 
ture of GRB hosts became somewhat more diverse; the offset of 
GRB-selected galaxies towards lower metallicities in the mass- 
metallicity relation ( [Levesque et al. 2010bp could, for exam¬ 
ple, be partially explained with the dependence of the metal¬ 
licity of star-forming galaxies on SFR (e.g., [Mannucci et al.j 
[2011 1 Kocevski & West|2011|l. Additionally, several metal-rich 


GRB hosts were discovered ( Levesque et al.|2010c[ Elliott et al. 
[2013[ [Schady et al.[[2015] l, and extensive observation in multi¬ 
band photometry revealed a population of red, high-mass, high- 
luminosity hosts, mostly associated with dust-extinguished af- 
terglo ws ([Kriihler et al. 201 1[ [Rossi et al.[[2012[ [Hjorth et al.l 
2012[ Perley et al. [2013b i. Similar GRBs were underrepresented 


in previous studies, illustrating the need for dust-independent. 
X-ray-selected samples such as The Optically Unbiased GRB 
Host survey (TOUGH, Hjorth et al.[2012|l or the Swift G RB Host 
Galaxy Legacy Survey (SHOALS, [Perley et al.|2015al l. 

Because of these inherent uncertainties from sample se¬ 
lection and small number statistics, the question of how di¬ 
rectly GRBs trace star formation, and how representatively they 
select star-forming galaxies remains a matter of debate (e.g., 
Michalowski et aL][2012[ [Vergani et al.[[20T5j [Schulze et al.[ 


2015[l: a metal-dependenc e on the GRB selection is shown by 
e.g., Graham & Fruchter P013[ l or Perley et al. ( [2013b[ l, while 


[Hunt et aTT ( [2014[ l advocate no strong evidence for GRBs to pro¬ 
vide a biased census of star formation. 

Elucidating the tight connection between star formation and 
GRBs also holds the promise of shedding more light on the na¬ 
ture of star-forming galaxies in general, especially at faint lumi- 
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Eig. 1; Brightness and redshift distribution of the galaxies in 
our sample. Magnitudes are primarily measured in the R/r Al¬ 
ter (94% of the sample). Eor the remaining hosts, V or ///-band 
magnitudes are shown. Upper limits come from events for which 
the spectral continuum was dominated by a bright afterglow, they 
are denoted by downward triangles and are not included in the 
brightness histogram. Galaxies that hosted dusty GRBs (sight¬ 
line A^^® > 1 mag) are indicated by red stars and red-filled his¬ 
tograms. 


nosities as well as their evolution at high redshift. Detailed spec¬ 
troscopic observation of GRB hosts, however, remained chal¬ 
lenging, in particular at z > 1; prominent tracers of the physical 
conditions in the hot gas are redshifted into the NIR where spec¬ 
troscopy traditionally is much less efficient. Spectroscopic data 
for z > 1 GRB hosts from emission lines is therefore available 
for only a handful of cases (e.g., Chen[[2012 Eriis et ar|[2015[ 
[Piranomonte et al.|2015[), and even at z < 1 the largest samples 


([Savaglio et al.[2009|[Levesque et al.[2010a[[Graham & Pruchter[ 


2013| l contain only 10-15 events with detailed information on the 
host’s gas properties. 

Here, we present initial results from emission-line spec¬ 
troscopy of 96 GRB-selected galaxies in the redshift range 0.1 < 
z < 3.6, with the bulk of the targets (76%) between 0.5 <z < 2.5. 
In this work, we focus on strong recombination and nebular lines 
such as the Balmer series and the forbidden transitions of [O ii] 
or [Om]. The Lya emission properties of the sample, a stack¬ 
ing analysis, and photometric follow-up observations will be dis¬ 
cussed in detail elsewhere. 

Throughout the paper, we report line fluxes in units of 
10“*^ ergcm“^ s“', magnitudes in the AB system, wavelengths 
and redshifts in vacuum and a heliocentric reference frame, and 
errors at Icr confidence levels. We assume concordance cosmol¬ 
ogy (Planck Collaboration [2014 = 0.315, Qa = 0.685, 

Hq - 67.3 km s“‘Mpc^Q, and a solar oxygen abundance of 
12 H- log(0/H) = 8.69 ( [Asplund et al^IOO^ . 


2. Sample and observations 

2.1. Targets 


To assemble our target list, we rely heavily on previous cat¬ 
alogs, in particular TOUGH, BAT6 ([Salvaterra et al. 2012) , 
from GROND ( [Greiner et al.[[20TTl [Kriihler et al.[[201 J), and 
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SHOALS for which we retrieve X-Shooteip]( |Vernet et al.]20TT 


optical/NIR spectroscopy from the ESO archive (Pis: de Ugarte 
Postigo, Flores, Fynbo, Kaper, Kriihler, Malesani, Piranomonte, 
Rossi, Schady, Schulze). In addition, we also use GRB afterglow 
observations (PI: Fynbo) obtained with X-Shooter in target-of- 
opportunity mode. 

We immediately removed events from the analysis that 
belong to the class of short GRBs such as GRB 130603B 
(e.g., |de Ugarte Postigo et al.|[2014[ ). We include the host of 
GRB 100816A even though it cannot be uniquely attributed to 
either long or short category ( [Oates et'ar|2010| l. To maintain a 
clean sample selected through genuinely long GRBs, this bor¬ 
derline case is not used in the discussion of galaxy properties in 
Sect.|5l GRB 10031 6D at z = 0.0592 (e.g., [Starling et al.|2011 


[Bufano et al.[[201^ is the closest GRB in the sample. As the 
X-Shooter slit covers only a small fraction of the host light, we 
cannot infer galaxy-integrated properties. We therefore also omit 
the GRB 100316D host galaxy in Sect.|^ 

After an initial reduction and screening of the afterglow or 
host spectra, we excluded those events for which neither emis¬ 
sion lines nor the stellar continuum of the host is detected: in 
the absence of a robust systemic redshift or in the presence of 
a very bright afterglow, no meaningful limits on emission-line 
fluxes can be established. We removed a total of nine host spec¬ 
tra from the sample in this step, primarily targeting very faint 
galaxies with R > 25.5 mag (GRBs 050406,060923A, 060923C, 
060919, 090926A, 101219B, 110709B). Two excluded spectra 
are likely those of foreground objects (i.e., misidentifications, 
GRBs 070808, 081210). Out of these nine host spectra, two 
have accurate afterglow redshifts (GRBs 090926A at z = 2.11, 
101219B at z ^ 0.55). GRB 060923A has a photo-z of z ~ 2.6 
( Tanvir et al.|2008| Perley et al.|2013b I. 

The total number of GRB host galaxy spectra presented here 
is 96, and nebular lines are detected for 91 of them. This sam¬ 
ple advances in two crucial aspects from previous works: For the 
first time, it extends significantly beyond redshift of order unity 
owing to the sensitivity of X-Shooter up to 2.5 fim. Our observa¬ 
tions thus cover prominent tracers (Ha and [N ii]) of the ionized 
gas up to z ~ 2.5, and [On], HjS, and [Oiii] up to z ~ 3.5. Sec¬ 
ondly, number statistics are approximately a factor five larger 
than in previous spectroscopic samples. A log summarizing the 
X-Shooter observations is provided in Table 
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Fig. 2: Comparison of the cumulative brightness distributions of 
the host of our X-Shooter emission-line spectroscopy (XES) in 
blue against representative GRB host samples in black and olive 
lines ( [Hjorth et al.[[2012[ [Vergani et al.[[2015) . In the top left 
of each panel, we show the K.S.-test p-value to reject the null- 
hypothesis that both samples are drawn from a similar parent 
distribution. 


(Fig.[T). It is hence shifted significantly towards lower redshifts 
when compared against GRB (afterglow) redshift distributions 


([Greiner et al. 201 f||Jakobsson et al.[2012[[Salvaterra et al.[2012[ 


[Perley et^ 2015a i, an obvious result of a selection bias favoring 

z < 3.5 galaxies for host spectroscopy. The brightness distribu¬ 
tion of the galaxies peaks at around R ~ 25 mag with a median 
value of Rtned ~ 24.3 mag. The basic properties of our galaxy 
sample, optical brightness and redshift are shown in Fig.[2 


2.3. Selection Effects 


2.2. The sample 


The presented X-Shooter observations were selected from the 
ESO archive, obtained by several groups and programs over the 
course of several years (2009-2015), and are hence diverse. They 
can however be divided into two main sub-categories: 

First, dedicated host spectroscopy, which aims on either 


measuring redshifts of GRBs for sample studies (e.g., Kriihler 
et al.|2012b [Salvaterra et al.|2012[l or the p hysical properties of 
individual GRB hosts (e.g., [Vergani et al.[[2011[ [IGiihler et al. 
2012a[l. Second, we also discuss GRB afterglow observations in 


which host emission lines are detected above the afterglow con¬ 
tinuum. 

Our GRB hosts are faint (21 mag < R < 26 mag with only a 
few exceptions) and their redshift distribution is broad extending 
up to z ~ 3.5. It has a median of Zmed = 1-56 and 35% of the tar¬ 
gets are at 0.5 <z< 1.5, and another 41% between 1.5 < z < 2.5 


* X-Shooter is a medium-resolution, cross-dispersed echelle spec¬ 
trograph mounted at ESO’s Very Large Telescope and sensitive from 
the atmospheric UV cutoff (~ 3000 A) up to the X-band. 


Galaxy samples obtained from GRBs are different than those 
from flux-limited surveys. Since the galaxy is localized by an 
external high-energy trigger - the GRB and its afterglow - the 
selection is initially independent of the galaxies dust content, 
brightness or SFR. If we studied the hosts of all (or a represen¬ 
tative subset of) GRBs, we could examine their properties irre¬ 
spective of any flux limit. 


As we, however, heterogeneously draw our GRBs from a 
larger parent sample of Swift GRBs, it is important to test how 
representative our GRB hosts really are. This implies under¬ 
standing the selection effects on the GRBs itself as well as those 
introduced through the follow-up observations. Dedicated spec¬ 
troscopic host observations favor, for example, optically brighter 
galaxies on average: a precise slit alignment requires a prior de¬ 
tection via broad-band imaging (Sect. [23(2]). 


A further critical parameter is the amou nt of visual attenua- 
tion by dust along the GRB sight-line ( Kriihler et al. 2011 
[Hjorth et al. [2012] [Perley et al.|2013b|[20l5't3l l because it corre- 
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lates with host properties - galaxies hosting dust}0GRBs have 
an order of magnitude higher stellar mass, luminosity, and SFR. 

The presence of GRBs with suppressed optical afterglows 
(e.g., |Groot et al.|1998[ [Levan et al.|2006| l constituted a signif¬ 
icant uncertainty in previous studies because they were under¬ 


represented in redshift distributions (e.g.. Fig. 9 in Hjorth et al. 


\2Q\2) . As X-Shooter is extremely efficient in providing galaxy 
redshifts ( Krithler et al.|2012b| l, there is an above average frac¬ 
tion of dusty (Ay*^" > 1 mag) GRBs in this work: it is 35 + 5 % 
here, but 20 % to 30 % in unbiased GRB samples ([Greiner et ah] 


2011[[Melandri et al.|2012[ Littlejohns et al.|2015| l in a similar 

redshift range (Sect. |2.3.3| l. 


2.3.1. Absolute statistics 


First, we examine absolute number statistics at z < 1, where the 
sample of GRBs with redshifts is highly complet^ In total, we 
have targeted ~50% (32) of all z < 1 GRBs in the Swift-era.. 
The remaining 50% were either too northern, too close to the 
Galactic plane, or not observable because of technical or tele¬ 
scope time limitations. From the initially targeted 32 galaxies at 
z < 1, we detect emission lines for 30 (94 %). The only two 
non-detections are GRB 110715A (z = 0.82), which spectrum 
was dominated by a R = 18.5 mag afterglow (jPiranomonte et al. 
[MTT) , and GRB 101219B at z = 0.55 ( |Sparre et al.|20TTF 


2.3.2. Brightness distributions 

To test whether the pre-imaging introduces biases towards opti¬ 
cally bright hosts, we compare the distribution of R-band mag¬ 
nitudes (Fig. with those of two re presentative and re dshift- 
complete samples ( [Hjorth et al.|2012[|Vergani et al.|2015] l which 
we use as a control group. 

Figure shows the respective cumulative brightness distri¬ 
butions in four redshifts intervals. At low redshift (z < 2), a 
K.S.-test returns no strong evidence that both samples draw from 
different parent distributions (K.S.-test p-values of 0.87 and 0.50 
at z < 1 and 1 < z < 2, respectively) - except possibly a ~20% 
fraction of very faint hosts in TOUGH at 1 < z < 2 that is 
not present here. Above z ~ 2, however, a selection effect is 
clearly evident: our galaxies are substantially brighter than the 
reference sample (K.S.-test p-values of p - 0.96 and p - 0.996 
at 2 < z < 3 and 3 < z < 3.8, respectively). 


2.3.3. The Fraction of Dusty GRBs 

Our z < 2 galaxies therefore have a similar brightness distribu¬ 
tion to unbiased GRB host samples. The over-proportionality of 


^We prefer to work with dusty GRBs instead of dark GRBs. The 
definition of dark GRBs is not unique ( jjakobsson et al.|2004[ [van derj 
[Horst et al.|2009| l, depends on the physics of the GRB shockwave and on 
the time and filter when the observations are performed ( [Greiner et al.[ 
[2011[ |. The rest-frame dust column density along the GRB sight line 
Ay*'” is thus a more physical quantity than the darkness of a afterglow. 

^Because of the exquisite X-ray positions from Swift and exten¬ 
sive follow-up efforts from the ground with e.g., TOUGH, GROND, 
or SHOALS, low-redshift GRBs are very unlikely to be missed if they 
have favorable observing conditions. 

‘^Representative, or unbiased, samples fully exploit the advantages 
of the high-energy selection by GRBs. Since they are based on criteria 
that are not connected to the burst’s physical environment, they do not 
depend on dust obscuration or galaxy luminosity. They allow us to study 
afterglows or hosts in a truly representative manner. 
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Fig. 3: Redshift dependent fractions with errors of dusty after¬ 
glows from our X-Shooter emission-line spectroscopy (XES) in 


blue against representative GRB samples in black (Greiner et al. 
20TT|[Covino et al.|2013[l. 


galaxies with dust-reddened afterglows in the sample, however, 
could nevertheless have an impact on their physical parameters. 

As representative comparison sample of the dust content to¬ 
wards GRBs ( A^*^^), we avail of the union of the afterglows from 
Greiner et al. ([201 1[) and Covino et al. (2013[l. After excluding 


duplicates, and limiting the samples to the same redshifts as the 
host spectroscopy (z < 3.6), 75 GRBs remain to constrain the 
fraction of dusty (A^^® > 1 mag) GRBs. Figure compares 
these values with those of our GRB0 

The ratio of dusty-to-normal afterglows is somewhat higher 
at z ~ 2 than at z < 1 (Figure]^. In addition, our sample has a 
slightly higher fraction of dusty afterglows at all redshifts than 
uniform GRB control samples. The data from galaxies selected 
through dusty GRBs are highlighted by stars in the following 
plots to illustrate their influence on correlations and GRB host 
evolution. 


2.4. A Stochastic Analysis 

To correct for the over-proportionality of dusty GRBs, we sta¬ 
tistically create samples that have the right fraction of dusty-to- 
normal afterglows as follows: We use the initial 96 hosts and 
randomly draw 10® sub-samples while using the previously es¬ 
tablished dusty GRB fraction as a strict prior. Specifically, this 
means applying a redshift-dependent probability to prevent a 
galaxy which hosted a dusty GRB to enter in sub-samples. 

The remaining analysis in this paper is performed with all 
10® sub-samples where each is a random representation of the 
input catalog but with a correct proportionality between dusty 
and normal GRBs. Nebular line fluxes are described by a Gaus¬ 
sian probability distribution with mean and standard deviation 
given by the raw measurement and error. Since the likelihood 
of keeping dusty GRBs has an error itself, we propagate in this 
way not only the statistical uncertainty due to sample size and 
the measurement error, but implicitly also the uncertainty in the 
prior, i.e., the dusty GRB fraction (28 + 8 % at 2 < z < 3, for 
example). 


^ Taken from the literature 

Cenko et aL|2009l|KrUhler et al.|201ll 

[Zafar et al.|2011| Perley et al. 

2015a[l or derived from our own after- 


constrain A™® > 1 mag or when the optical/NIR afterglow is signif¬ 
icantly underluminous with respect to the X-ray data. For example, 
GRB 120119A is dusty, although it had a bright optical counterpart 
([Morgan et al.[[2014[>, but GRB 100606A is normal (not-dusty), even 


though it did not jNicuesa et al. 2010 1 . 
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Fig. 4; Probability distribution of chance coincidence, or the 
likelihood of having a given number of galaxies in the total sam¬ 
ple wrongly associated with the GRB. The probability that we 
have misidentified a total of 2 out of 96 galaxies is, for example, 
30%. 


The result of this procedure are 10^ different distributions of 
host parameters, correlations, and best-fit values. In the follow¬ 
ing, we then provide for all parameters or correlations the me¬ 
dian of the a-posteriori probability and the Icr-equivalent range 
as error to represent their distribution. 

Whenever we study GRB host parameters, we will plot two 
data sets: First, the more meaningful distribution restricted to 
z < 2 (Sect. |2.3.2| l and derived from the stochastic sample in¬ 
cluding error-bars always in blue. Second, we will also show the 
histogram of all measured values in black to illustrate the impact 
of the sample cuts. At z > 2, we show the resulting parameters 
and distributions (without error regions) as reference. The selec¬ 
tion effect towards optically bright galaxies at z > 2, however, 
will often prevent us from deriving robust conclusions in this 
redshift range. 


brightnesses f rom Fig. [T| to assess how many similarly bright 
field galaxies ( Capak et al.||2004] [Gawiser et al.||2006| l are sta¬ 
tistically expected within an area corresponding to the GRB’s 
localization accuracy. After excluding those events for which 
the afterglow fine-structure/host redshift match, we calculate the 
probabilities for a given number of misidentifications through a 
Monte-Carlo method with 10® trials. The resulting distribution is 
shown in Fig. and has a peak at around 1-2 events, with a tail 
out to a total of 5 misidentified galaxies. Very few field galaxies 
are thus interloping in our GRB host sample and their impact on 
the main conclusions will be limited.X-Shooter 


3. X-Shooter Optical/NIR spectroscopy 

X-Shooter operates in three dichroic-separated arms simultane¬ 
ously: the ultra-violet, blue (UVB, 3000 A-5600 A), the visual 
(VIS, 5500 A-10 020 A), and the near-infrared (MR, lOOOOA - 
24 800 A) arm. Each of the three arms has its own slit, echelle 
grating, cross-disperser, and detector, and we refer to |Vernet| 
|et al.| ( |201 1[ ) for a detailed description of the instrument. The 
resolving power of X-Shooter depends on the arm and slit width 
and is^pically between R - 4000 and R - 10000 for the used 
setup^ Unless dictated by other observational constraints such 
as nearby bright objects, the observations were obtained in a 
nodding pattern in multiples of an ABBA sequence. The nod- 
length is constrained by the short slit-length of X-Shooter (11'.'O) 
and is usually 5'.'0. 

Whenever possible, slit loss due to atmospheric dispersion 
was minimized by observing at the parallactic angle and at min¬ 
imum airmass. A large subset of the host observations were ob¬ 
tained before X-Shooter’s ADC failure in August 2012 and most 
of the emission lines are located in the NIR part of the spectrum 
(e.g., the bluest of the lines of interest ([On]) is already above 
10000 A at z > 1.5) where atmospheric dispersion is less severe. 
The lack of ADCs was not a primary concern in the analysis. 


2.5. Misidentification of Hosts 


3.1. Data reduction 


We use standard procedures ( [Bloom et al.||200^ jPerley et al.j 
[2009 j l to assign galaxies to GRBs by using the afterglow position 
with the highest accuracy available. A crucial basis of this work 
are accurate GRB positions which are routinely available only 
since the launch of Swift (jGehrels et al.||2004l [Burrows et al.j 
[2005^ . 

In the best cases (16 of 96 GRBs), an afterglow redshift is 
set by fine-structure lines (e.g., [Vreeswijk et al.[[2007| l, and a 
matching emission line redshift as well as a sub-arcsecond error- 
circle leave hardly any room for ambiguity. In total, 73 of the 
GRBs have localizations with an accuracy of around (y.'5 through 
radio, optical or Chandra X-ray afterglow observations. 

In the remaining cases (23 of 96 GRBs), the host association 
relies on the X-ray position from Swift which constrains the GRB 
position to a circle with ~1'.'5 radius at 90% confidence. The 
probability of misidentification, i.e., the galaxy that we study is 
not related to the GRB of interest is then a function of positional 
accuracy, galaxy brightness, and density of galaxies for the given 
galaxy type, brightness, and coordinates (e.g., [Rossi et al.[2014|l. 
For average values of galaxy densities in enipty fields |Capa^ 
et al.[[20()^ Gawiser et al. 2006[), positional errors of T.'5, and 


galaxy brightnesses of r ~ 25th magnitude, the chance of finding 
a random galaxy in the error circle is around 1-2%. 

To quantify the total number of chance coincidences, or 
misidentifications in the full sample, we use the GRB host 


The basis of the data reduction was the X-Shooter pipeline sup¬ 
plied by ESO in its version 2.2.0 or higher ( [Goldoni et al.[2006[ 
[Modigliani et al.[|2010| l, which we used for flat-fielding, order 
tracing, rectification, and initial wavelength calibration and flex¬ 
ure compensation via arc-lamp frames (also Sect. [3.2[ ). These 
steps were applied to the individual frames of the individual arms 
of each set of spectra. During rectification, we chose a dispersion 
of 0.4 A/px (UVB/VIS arm) and 0.6 A/px (NIR arm). This min¬ 
imizes correlated noise, while at the same time maintaining suf¬ 
ficient sampling for emission lines down to a velocity dispersion 
cr of ~ 20kms“'. 

We use our own software and algorithms for bad-pixel and 
cosmic ray detection, as well as sky-subtraction and stacking of 
individual exposures. We further corrected the inter-order back¬ 
ground for an apparent problem in the X-Shooter pipeline. In the 
UVB and VIS arm, we estimate the sky background locally in a 
small region around the spectral trace. Eor the NIR data with its 
abundant and strong sky-lines and high background, the frames 
for each set of observation are grouped in pairs of two such that 
they were taken in a different sky-position but as close in time as 
possible. Each two frames were then subtracted from each other 
for the background estimation. 

*’https://www.eso.org/sci/facilities/paranal/instruments/ 

xshooter/inst.html 
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A typical observation (Table [^l consisting of four expo¬ 
sures in the UVB/VIS, and twelve in the NIR arm leads to four 
and six sky-subtracted, rectihed, wavelength and flux-calibrated 
(Sect. |3.3[ ) frames in the UVB/VIS and NIR arm, respectively. 
These single frames are then averaged using a weight function 
from the variance or the signal-to-noise ratio (S/N). Previously 
detected bad-pixels and cosmic rays are masked in the stacking 
process. 

From the combined two-dimensional frame we then opti¬ 
mally extract the one-dimensional spectrum. The weight func¬ 
tion of the extraction is derived using the collapsed profile of 
either the emission lines themselves or the spectral continuum in 
case it is detected at a sufficient S/N level. 


3.2. Wavelength calibration 

X-Shooter is mounted on the VLT Cassegrain focus, so it suf¬ 
fers from instrument flexures and misalignments of the wave¬ 
length solution between day-time calibration and night-time sci¬ 
ence observations. We correct for these flexures first by using the 
ESO pipeline and frames supplied by X-Shooter’s active flex¬ 
ure compensation. In a second step, we use a theoretical emis¬ 
sion spectrum of the Paranal sky ( |Noll et al.|2012| l. In the VIS 
and NIR arm, we divide the spectrum into 10-15 slices, and 
cross-correlate the observed and theoretical sky-spectrum. This 
provides wavelength offsets averaged over the respective wave¬ 
length interval. A low-order polynomial was then fitted to the 
data, and the respective wavelength offset as derived from the fit 
was applied to the wavelength scale. Typical offsets are smaller 
than around 0.4 A in the NIR arm and < 0.1 A in the UVB/VIS. 
Because of the lack of strong sky features in the UVB arm blue- 
ward of 5000 A, we apply only a constant offset to the wave¬ 
length scale in the UVB arm. 


3.3. Flux calibration 

Correct flux scales of the emission line spectra are fundamen¬ 
tal because they have a strong impact on the measurement of 
global galaxy characteristics. We first flux-calibrated each of the 
individual spectra with the nightly standard star taken with a slit- 
width of 5'.'0. Standard stars are different for each observation, 
and the full set is given at the ESO webpage^ The resulting 
response curves were thoroughly validated against each other. 
To limit the effect of a time-dependent illumination by a time- 
varying flat-field lamp, particularly evident in the red part of the 
NIR arm, we used the same flat field for science and standard 
frames to correct for pixel-to-pixel variations. 

Even with a well-defined instrumental response, the finite slit 
width inevitably leads to slit-loss. Slit-loss depends primarily on 
the ratio between slit-width and the convolution of the point- 
spread function delivered by the telescope with the spatial extent 
of the sources on the sky. In addition, the sky transparency varied 
significantly during the course of some of the observations due 
to clouds, leading to an additional decrease in flux compared to 
the standard star observation. 

Accurate broadband photometry, however, is available to us 
providing means to estimate galaxy-integrated line-fluxes. The 
ratio between photometry and a synthetic brightness of the spec¬ 
trum integrated over the band-pass of the respective filter yields 
a scale factor to estimate slit-loss (which we call t). In the 
ideal case, multiple independent photometric bands in each of 

’http://www.eso.org/sci/facilities/paranal/instruments 

/xshooter/tools/specphot_list.html 
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X-Shooter’s arms allow us to corroborate and check the intra- 
and inter-arm flux calibration. Eor many hosts, this is not the 
case and we assume that t is constant in each arm. We then mul¬ 
tiply the line flux by t and propagate errors accordingly, t is in 
the range between 1 and 3, and its distribution has a median and 
Icrbounds of Tmed = 1 - 6 ^o 5 - 

Deriving t for the individual arms is only possible for a 
bright continuum (r < 23 mag), i.e., luminous hosts or after¬ 
glows. Eor most of the fainter hosts (r ~ 25 - 26 mag), the con¬ 
tinuum is detected in the bluer arms only leading to relatively 
large errors in t and the line fluxes. In cases where no contin¬ 
uum is detected, or no comparison photometry is available in the 
given arm, we assume that t is similar between all arms. In cases 
where we can measure t for different arms or wavelengths, it is 
consistent within errors. This provides confidence that the scal¬ 
ing procedure is reliable within the quoted error bars. 

Errors in t come from the photon noise of the spectral contin¬ 
uum, the accuracy of the comparison photometry, and the scatter 
of different filters in each arm if available. Line ratios of adja¬ 
cent lines (such as Ha and [N ii] or the [O iii] doublet and HjS) 
and the physical properties that they are proxies for, have the 
same scaling factor in all cases. We associate the line-flux error 
from photon statistics with a statistical error, and the one de¬ 
rived from the scaling factor to a systematic one, and provide 
full details of both errors in the line-flux tables (Tables]^ and [^. 
All fluxes, spectra, and photometry are corrected for the Galactic 
foreground reddening according to Schlafly & Einkbeiner| ( |201 1) 
assuming Ry - 3.08 ( |Pei|1992| ). 


3.4. Flux limits and sensitivities 


To better understand the fundamental properties and constraints 
of our observations, we look closer at the sensitivity limits. Erom 
a typical observational setup with an integration time of 4x600 s, 
we derive the noise characteristics of the spectrum. Artificial 
emission lines of Ha and the [O ii] doublet at various reds hifts 
with a typical intrinsic line width of cr = 70kms“' (Sect. 5.5 1 
are then added to the noise spectrum. 

We record the emission-line flux that would have been de¬ 
tected with a combined significance of at least 6 cr over the noise 
level in the one-dimensional spectrum. The significance level is 
set conservatively because we have no a priori knowledge of the 
redshift in many case^ Line-flux sensitivities for Ha and [O ii] 
are then converted into SERs using the methods from Sect. |4.4| 

Eigure [^illustrates the redshift dependence of the detectable 
SER. At z < 0.5, Ha lies in the visual arm of X-Shooter with 
its sensitive CCD and low sky-background. The SER limit is a 
few 10 “"^Moyr“'at z ~ 0.1 and increases to 8 ■ lO^^Moyr^'at 
z ~ 0.5. In the non-extinguished case (solid lines in Eig. |^, the 
SER sensitivity stays below 1 Mq yr“^ up to z ~ 1.7. 

Above z ~ 2.5, the sensitivity to unextinguished SER is 
driven by [On] and stays at around a few Mq yr^'up to z ~ 3.7 
above which it rises beyond lOMoyr^'. Adding dust increases 
these limits considerably, in particular for [O ii] (dashed lines in 
Eig.jg. 

Noteworthy are redshift ranges where the SE-tracers are lo¬ 
cated in regions of prominent telluric absorption bands (black 


*The look-elsewhere effect: It is for example very likely to find 3 cr 
noise peaks in the medium-resolution spectrum of X-Shooter just be¬ 
cause the parameter space to be searched is so large. If the redshift is 
known, the locations of the emission lines are fixed, and the probability 
of finding significant noise peaks by chance in the correct wavelength 
range drops substantially. 
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Fig. 5; Colored lines represent the sensitivity for Her and 
[O It] converted into a SFR. Solid/dashed lines denote the case of 
Ay - Omag and Ay - 2 mag, respectively. Black lines indicate 
redshift ranges where the line is located in windows of negligi¬ 
ble atmospheric transmission. Skylines result in high SFR limits 
on small scales, but have been omitted to enhance clarity. 


lines in the top panel of Fig. [^. In a small window around z ~ 1 
and redshift ranges of 1.7 < z < 2.0, 2.6 < z < 3.0, and z > 3.7, 
atmospheric absorption makes the detection of emission lines 
challenging from the ground. 


3 . 5 . Auxiliary photometry 


For most of the targets the photometry used to flux-calibrate the 
spectra is either provided by TOUGH ( [Hjorth et al.||2012| ), or 
the literature referenced in Table [T] If not available elsewhere, 
dedicated observations were performed with GROND (|Greiner| 
[erai:][2008ll at the 2.2 m M PG teles cope on LaSilla, HAWK - 
I (|Kissler-Patig et al.||20081l, FORS2 (|Appenzeller et al.||1998|l, 
LRIS at Keck joke et al.|l^ or ALFOSC at the 2.5 m Nordic 
Optical Telescope (NOT). 

The optical/NIR imaging data were reduced in a standard 
fashion using dedicated routines in pyraf/IRAF ( |Tody|[T99^ in 
the case of GROND, VLT, and NOT photometry or IDL for 
Keck/LRIS. Data reduction and photometry follow the proce¬ 
dures established in |Kruhler et al.|p008| l in the former case and 
is taken from the Keck GRB Host Survey ( |Perley et al.|2013^ i 
in the latter case. The photometric solution of each field was tied 
to the SDSS survey for those GRB fields which were covered 
by SDSS DR8 ( |Aihara et al.|2011) , zero-points derived from the 
observation of standard stars, or our own set of secondary cali¬ 
brators in the field. Typically, secondary calibrators were again 
derived from observing SDSS fields close in time to the GRB 
observations. 


4. Analysis and results 

4.1. Line flux measurements 

Fluxes of emission lines were measured by fitting a superposi¬ 
tion of Gaussian functions to the data, with the continuum set in 
small (+40 A) regions around the wavelength of interest free of 
emission-, telluric-, and sky-lines. In most cases, a single Gaus¬ 
sian and constant continuum describes the line with sufficient 
accuracy (Fig. [^. We also measured the line-flux through nu¬ 
merical integration, and its result was used in cases where the 
lines show clear velocity structure. 


For weak emission lines such as [Ne iii] or [N ii], we tied the 
Gaussian centroid (set by the galaxy redshift) and the line width 
to those of strong lines. 

Finally, a synthetic Paranal sky radiance and transmission 
spectrum ( [Noll et JI]|2012| l was used to flag regions of telluric 
absorption and sky emission, and affected data were automati¬ 
cally discarded from the line fit (Fig. |^. The results of the line 
flux measurements are given in Tables^and]^ 

Many of the emission lines are not detected at high signal-to- 
noise ratios (S/N). If the spectral range is covered by X-Shooter, 
we decided to quote the actual measurement instead of providing 
an upper limit to maximize the information on e.g., line ratios 
under the given photon statistics. 

We have not corrected the Balmer line fluxes for the under¬ 
lying Balmer stellar absorption even though we note that it can 
have a measurable influence on the line fluxes and thus the de¬ 
rived physical properties of the emitting gas ( Wierserna||20 11 | l . 
The S/N of the emission lines and stellar continuum is not high 
enough to resolve and fit the absorption component separately. 
Generally, the Balmer absorption correction is small for galaxies 
with masses below 10'° Mq. 

At z ~ 0.8, for example, it has a n equivalent width (EW) 


of 1 A for galaxies of IO'^Mq at H/I (iZahidetak 


2011 1 . The 


EW of the Ha line is high enough (EWrest > 20 A) to make 
the correction negligible, but it can have an effect on the other 
Balmer line-fluxes. We estimate <5% for H/3 on average, well 
below the average statistical error of the line flux measurements 
(Tableig. 

In cases where [Oiii]/I5007 is not well constrained by 
the data, but [Oiii]/14959 is, or vice versa, we assume that 
/[Oiii(.i 4959 )] = V3 •/[ 0111 ( 45007 )] ( |Storey & Zeippen|2000) l. 


4.2. Emission-line widths 

We measure the intrinsic line-width (cr) of the hot gas using neb¬ 
ular lines of highest S/N, usually Ha, or [O iii]. If the line-shape 
is adequately represented with a Gaussian function, we use the 
instrumental resolution AV and subtract it quadratically from the 

fitted EWHM, giving cr as cr = VfWHM^^^A^/( 2 V2ln2). 
If the line-shape is clearly asymmetric, we directly measure 
EWHM in the data. Errors on cr are derived via the error- 
spectrum and propagated accordingly. 

The physical line-width is well-resolved through the medium 
resolution of X-Shooter (typically AV ~ 35 km s“^ in the visual 
arm, or AV ~ 50 km s ' in the NIR arm). The line-width mea¬ 
surements are provided in Table 


4.3. New GRB redshifts 


To derive GRB redshifts from host galaxy data (Tabl e [5]l, w e 
use both the emission lines as detailed in Kriihler et al. (|2012bll, 
spectral breaks in the stellar continuum ( Jakobsson et al.|2012|l 
as well as standard photo-z techniques in LePhare (Amouts et al. 
|1999[ Ilbert et al.|2006| l from supplementary multi-band photom¬ 
etry. Eor galaxies with more than one detected emission line, the 
redshift solution is accurate, unique, and robust. Eor the other, 
somewhat more ambiguous cases we provide more information 
in Appendix [A| 

In particular, we provide redshifts for six GRBs 
(GRB 050714B at z = 2.438, GRB 06030^ at z = 1.56 


^The redshift of GRB 060306 has been discussed in ijakobsson et al.l 
\2Q12) and |Perley et al.| ( |2013b^ , already suggesting z = 1.55 as the 
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Fig. 6: Examples of line-flux measurements for nebular lines at different signal-to-noise level. In each panel, the background is 
the two-dimensional spectrum, where the spatial direction is plotted against the right y-axis. The color coding is indicated by the 
color bar. The one-dimensional spectrum with a Gaussian fit are shown in black points and line. Gray data points are excluded from 
the fit because of telluric or skyline contamination. A telluric transmission spectrum is plotted in a thin black line at the very top 
of each panel. From top left to bottom right: [On] (GRB 140506A, with afterglow background), [Neiii] (GRB 100424A), [Out] 
(GRB 070129), Hy6 (GRB 051022), Ha (GRB 070110), [Nii] (GRB 050416A). 


GRB 060805at z = 2.363, GRB 060923B at z = 1.509, 
GRB 070224 at z = 1.9 92, and GRB 0703 28 at z = 2.063) 
in the TOUGH sample ( [Hjorth et al. |2012[ ). This brings the 
redshift completeness of TOUGH to 87% (60 out of 69 events), 
with a median/average redshift of Zmed = 2.10 + 0.18 and 
Zavg = 2.20 ± 0.17, respectively. GRB 070328, GRB 070521 
(z = 2.087), GRB 080602 (z = 1.820), and GRB 090201 
z = 2.100 are part of BAT6, so its completeness is now 97% (56 
out of 58 events) with Zmed = 1.67 + 0.15 and Zavg = 1.75 + 0.13. 

In addition, we measure or constrain the redshifts of seven 
GRBs based on emission lines or spectral breaks in the stellar 
continuum (Table [^. 


4.4. Scaling relations for star formation rates 


Since Ha, the most reliable tracer of star formation in the 
rest-frame optical wavelength range, is not visible in all of 
the spectra, we use our sample to derive scaling relations be¬ 
tween the luminosity of other strong nebular lines and the dust- 
corrected SFRhc. To establish the comparison sample, we avail 
of those galaxies for which both Ha and Eb-v from the Balmer 
decrement (Sect. |5.3| l can be constrained. We convert the dust- 
corrected Ha luminosity into a SFR based on the widely used 
Kennicutt (19981 relation adopting a Chabrier ( 2003| l IMF. This 
provides us with a reference sample with reddening-corrected 
SFRs based on Ha. 

We then plot other strong emission lines against SFR (e.g.. 
Fig. 0 and minimize the scatter of the data against scaling re¬ 
lations in the form of log(5 URhu) - a + b ■ log (L 42 ). Here, 
lO'* is the SFR in Moyr ' for a benchmark line luminosity of 
L 42 = lO'^^ergs^', while b takes into account a possible evo- 


likely redshift based on a VLT/FORS spectrum in the former, and ex¬ 
tensive photometry in the latter case. The detection of multiple emission 
lines unambiguously confirms z = 1.56 as the correct redshift. 

*®The association of GRB 060805A is ambiguous because two 
galaxies are detected in the Swift XRT error-circle. Here, we provide 
the redshift for the fainter object. The brighter object has a redshift of 
z = 0.603 Oakobsson et al.|2012b. 


Article number, page 8 of 33 


Table 5: New GRB redshifts 


GRB host 

Redshift 

Spectral features 

GRB 050714B 

2.4383 

[Oiii](5007), Ha 

GRB 060306'= 

1.5585/ 1.5597 

[On], Ha 

GRB 060805A* 

2.3633 

Hj6, [Oiii](5007), Ha 

GRB 060923B 

1.5094 

[Oiii](5007), Ha 

GRB 070224 

1.9922 

[Oiii](4959,5007) 

GRB 070328 

2.0627 

[On], H/3, [Out] 

GRB 070521 

2.0865 

Ha 

GRB 080602 

1.8204 

[On], [Out], Ha 

GRB 090201 

2.1000 

[Oii],HyS, [Out], Ha 

GRB 100508A 

0.5201 

[Oii],HyS, [Out], Ha 

GRB 111129A 

1.0796 

[On] 

GRB 12021 lA 

2.4+ 0.1 

Fy-a break 

GRB 120224A 

1.1 +0.2 

Balmer break 

GRB 120805A 

3.1 +0.1 

Fy-a break 

GRB 121209A 

2.1 +0.3 

Balmer break 

GRB 140114A 

3.0+ 0.1 

Fy-a break 


Notes. The horizontal lines denote the separation between TOUGH 
GRBs (upper part), BAT6 (middle part) and others. Emission lines 
of the host of GRB 060306 have a double-peaked profile, and we thus 
provide two redshifts. Ambiguous association because two galaxies 
are present in the XRT error-circle. This is the redshift for the fainter of 
the two objects, for details see|Jakobsson et aL](|2012|l. 


lution between SFR and line flux with galaxy properties. These 
scaling relations then constrain the SFR for hosts for which other 
nebular lines than Ha are available. 

A slope of b - 1 provides a direct proportionality between 
e.g., [O ii] luminosity and SFR ( |Kennicutt||1998[ |Savaglio et al.| 
2009|l. b significantly different from unity indicates that the rela¬ 


tion between SFR and line flux changes with galaxy properties, 
for instance metallicities, dust attenuation or ionization (e.g., 
Kewley et al.| 20041. One would thus ideally add a further de¬ 


pendence in the relation based on e.g., metallicity, stellar mass 
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Fig. 7; Scaling relations between emission-line luminosity and 
SFR. Uncorrected/de-reddened [O ii] luminosity against SFR in 
the upper and lower panel, respectively. Error-bars in the lower 
panel are omitted to enhance clarity (typical errors in the top 
left corner). The color-coding of the individual measurements 
corresponds to Ay as indicated by the color bar. The solid black 
lines show the best-fit scaling relations from Equation The 
gray-shaded area denotes the Icr scatter in the relation - 0.45 dex 
in the upper and 0.19 dex in the lower panel. 


with a dispersion of 0.25 dex. Eor Hj8 the equivalent relation is 

log(5 FRna) = 1.47 + 0.98 ■ log (L(H/3,™eor,42)) (3) 

with a scatter of 0.36 dex from 57 galaxies. Eor the Balmer lines, 
dust-corrected relations are by construction equivalent to the one 
of |Kennicutt| ( | 1 998] l. 

Eor the [O iii](/l5007) line, the scaling relations are 

logiSFRua) = 0.98 -t 0.81 ■ log(L[OlIl],uncor,42) 

log(5 FRna) = 0.71 -H 0.92 ■ log (L[0 IIl],dered,42) (4) 

with a dispersion of 0.42 dex for the uncorrected, and 0.28 dex in 
the dereddened case (from 57 individual galaxies in total). The 
relatively large dispersion, even after removing the effect of dust 
reddening, illustrates that [Out] is the least reliable SE-tracer 
considered in this study. In the following, we use the emission¬ 
line that yields the smallest uncertainty in SER when propagat¬ 
ing Icr dispersions accordingly when necessary (typically first 
Ha, then Hj3, [On], and [Oin] only if no other emission line is 
detected). 

It is immediately clear that the scaling relations detailed 
above depend on the distribution of galaxy properties in the ref¬ 
erence sample. Or, in other words, we consider these scaling 
relations most useful and accurate for low-mass, high-redshift, 
star-forming galaxies which have the average color excess of our 
sample (Ay ~ 0.8 mag or ~ 0.25 mag. Sect. ig. 


or luminosity for accurate SERs (Kewley et al. 2004 


Gilbank 


|et al.|2010| ). In the absence of further information on the galaxy 
apart from line flux, however, all the physical changes in the SER 
to line-flux relation are implicit in the empirical factor b. 

[O It] for example, is strongly correlated with SER (Eigure|7]i 
and based on 60 hosts at 0.3 < z < 2.5 we derive 


log(SFRHa) = 1.02 -H 1.07 ■ log(L[Oll],uncor,42) 
log(5 FRna) = 0.63 4- 1.03 ■ log (L[0 Il],dered,42) , (1) 


where T[Oii],uncor,42, T[Oii],dered,42 are the measured and dered¬ 
dened [O It] luminosity in units of 10"^^ erg s ', respectively. The 
Icr dispersions in these relations are 0.44 dex for the dust- 
uncorrected and 0.19 dex for the dust-corrected luminosity. Even 
though the scatter is much lower in the latter case, the former is 
useful in most practical cases. Eor a reliable dust-correction, the 
measurement of at least two Balmer lines is required, which in 
itself provides a much more reliable SER than through Equa¬ 
tion [1] 

Eigure also illustrates the limitation of this procedure. 
Since [O n](/13727 A) is strongly attenuated by dust, there is con¬ 
siderable scatter in the relation using the dust-uncorrected lu¬ 
minosity. Even though the calibration sample contains galax¬ 
ies over nearly full range of metallicity probed by the obser¬ 
vations (8.2 <12-1- log(0/H) < 8.9), the uncorrected Equa¬ 
tion will significantly under-predict the SER for metal-rich 
(12 4- log(0/H) > 8.7) galaxies because of their high dust con¬ 
tent. If a reliable dust-correction is available, the dust-corrected 
equation [T] then provides SERs with an average dispersion of 
0.2 dex over the full metallicity range probed. 

Eor those cases where no dust correction can be performed, 
i.e., only one Balmer line is detected in the spectrum, we de¬ 
rive for Ha (based on 61 individual galaxies with 8.2 < 12 4- 
log(0/H) < 8.9) 

log(5 FRna) = 0.92 4- 1.04 ■ log (T(Ha,uncor, 42 )) (2) 


5. Physical properties of the GRB host population 

5.1. GRB hosts in the BPT diagram 

Eigure|^shows the GRB hosts in the Baldwin-Phillips-Terlevich 
(BPT) diagram ( [Baldwin et al.|[T98T] ) that allows emission-line 
objects to be separated into star-forming galaxies and AGNs 
based on the [Nii]/Ha and [Oiii]/i5007/Hy8ratios (e.g., |Kewley| 
et aLpOOT). In our sample, 34 GRB hosts have adequate con¬ 


straints on the four emission-line fluxes to establish their posi 
tion in the BPT diagram. As we are interested in line-flux ratios, 
we use here (and in the following whenever line-ratios are im¬ 
portant) only the statistical errors from Tables |^and[^ 

The redshift range of the GRB hosts extends from z = 0.3 
to z = 2.5. As background sample, we plot the logarithmic 
space density of 330000 local (z ~ 0) SDSS (jAbazajian et al. 


et al. 


2009 1 galaxies from the MPA-JHU catalogue (e.g., Kauffmann 
2U03] l, 11 as the redshift-dependent classification lines 


between AGN and star-forming galaxies as gray dashed lines 
([Kewley et al.|2013b) l. The ridge line of SDSS galaxies ( jBrinch-j 
|mann et al.|2008| l isshown in black. 

Our GRB hosts occupy a different phase-space than SDSS 
galaxies in the BPT diagram because they are mostly above the 
ridge line which denotes the highest density of local star-forming 
galaxies. A similar offset is well-observed in many high red- 
shift galaxy samples (e.g., |Shapley et al.|2005[|Kriek et al.|2007[ 
Steidel et al.|2014)l and also for galaxies hosting H-poor super- 


luminous supernovae ( jLeloudas et al.|20T5l ). This offset is often 
attributed to harder ionization fields, higher ionization parame¬ 


ters or changes in the ISM properties (e.g., Brinchmann et al. 


|2008l|Kewley et al.|2013al|Steidel et al.|2014|. 

There is some evidence of an evolution of the position of 
the targets in the BPT diagram with redshift. The top and mid¬ 
dle panel of Pig. show the minimum difference to the SDSS- 
based ridge line (A([Oiii]/Hj0)) as well as [Oiii]/H/3. On aver¬ 
age A([Oiii]/Hj0) and [Oiii]/H;6 increase by 0.11 + 0.02 and 
0.08 + 0.03 dex per unit redshift. These values are derived as 
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Fig. 8; Baldwin-Phillips-Terlevich (BPT) classification diagram 
for GRB hosts. The background sample is the logarithmic space 
density of local (z ~ 0) SDSS galaxies with line-fluxes from 
the MPA/JHU catalo^^ Colored circles are 34 individual GRB- 
selected galaxies with redshifts indicated by the color bar. The 
increasingly darker dashed gray lines are the redshift-dependent 
differentiation lines between star-forming galaxies and AGN 
of Kewley et al. ( |2013b| l at the indicated redshifts. The black 
solid line is the ridge line, i.e., the line with the highest den¬ 
sity of SDSS galaxies. Dusty/normal GRBs are indicated by 
stars/circles. 
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Fig. 9: Minimum difference of the GRB hosts from the ridge-line 
(illustrated by the dots) of z ~ 0 SDSS star-forming galaxies (top 
panel). We show [Oiii]/Hy6 in black and ([Oiii]-i-[Oii])/Hy6 in red 
as a function of z in the middle and bottom panel, respectively. 
In all panels, solid lines denote the most-likely linear regression, 
with regions of 68 % probability in shaded areas. Dusty/normal 
GRBs are indicated by stars/circles. Error bars are omitted to 
enhance clarity and replaced by typical errors in the bottom-right 
corner. 


the median of all stochastic samples (Sect. |2.4] l with errors repre¬ 
senting the region of 68 % likelihood of all trials. The fractions of 
samples producing no evolution in [Oiii]/H/3, and A([Oni]/HjS) 
are p = 0.015 and p - 0.0018, respectively. A rising trend is 
not seen (bottom panel of Figure 0 in the ratio of the summed 
oxygen line flux ([Out] H-[Oii])to Hj8 (0.01 + 0.04 dex per unit 
redshift). 


Changes in the metallicity of the nebular gas are unlikely 
to play the dominant role in the displacement between high- 
redshift and local galaxies as this would move galaxies primarily 
along the ridge line. Photoionization models from [Steidel et al.| 
( |2014| l, for example, predict very little contribution from metal¬ 
licity changes to the observed offset in the BPT diagram. In addi¬ 
tion, the metallicity distributions of the sample are comparable at 
z < 1 and 1 < z < 2 (Sect. [5^6^ , indicating that the larger differ¬ 
ence of galaxies at higher redshift is not driven by a decreasing 
gas-phase oxygen abundance. 


Our observations are thus consistent with the hypothesis that 
the difference in the location in the BPT diagram between GRB 
hosts and z ~ 0 star-forming galaxies is caused by an increase 
in the ionization fraction, i.e., for a given metallicity a larger 
percentage of the total oxygen abundance is present in higher 
ionization states at higher redshifts. This could be caused by a 
harder ionization field originating for example from hot O-type 
stars ( |Steidel et al.||2014] l that emit a large number of photons 
capable of ionizing oxygen into [Out]. 
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5.2. Ionization evolution 


The ratio of [Oiii]/[Oii], which we plot in Fig. 10 against red¬ 
shift, Ay, and SFR, is often used as a proxy for the ionization pa¬ 
rameter. A rising [O iii]/[0 ii] ratio with redshift is evident from 
79 events with adequate wavelength coverage to constrain both 
[Oiii] and [On]. [Oiii]/[Oii] increases from 0.9 at z ~ 0.3 by a 
factor of 3 to ~ 3 at z ~ 3.3. The fraction of iterations return¬ 
ing no correlation between the [O in]/[0 ii] ratio and redshift is 
p = 0.0018 for the full sample and p = 0.0039 when restricting 
the redshift range to z < 2.7. 

The predicted [O ni]/[0ii] ratio at z = 3.3 is remarkably sim¬ 
ilar to the ratio of the composite LBG spectrum at z ~ 3.5 
( |Maiolino et al.|[2008[ [Troncoso et al.|[T014| l. The evolution in 
[O in]/[0 It] is again unlikely to be caused by metallicity effects 
(lower metallicity galaxies have higher [O iii]/[0 ii] ratios) be¬ 
cause there is no statistically significant difference seen in metal¬ 
licity between the z < 1 and 1 < z < 2 sample (Sect. ig. 

Even though both the ionization parameter and SER are cor¬ 
related with redshift, the bottom panel in Eig. 10 shows no ob¬ 
vious relation between [Oiii]/[Oii] and SER (Sect. I5.4|l. It does 


show (middle panel in Eig.[T0|) an anti-correlation with dust con¬ 
tent (Av = 0.35 - 0.90 ■ logRO III] /[On]), with p = 0.037). 
The dust contenj^as parametrized through a color-excess Eb-v 
or visual attenuation Ay of a galaxy (Sect. correlates with 


* * http://www.mpa-garching.mpg.de/SDSS 

'^Here, and in the following, Eb-v or Ay always refers to the rest- 
frame color excess/visual attenuation at the redshift of the target galaxy. 
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2.0 ■ 


1.0 



analysis (Sect. [233] and Sect. |2.4| l. The black histogram repre¬ 
sents all raw measurements without corrections. The lower panel 
shows the cumulative distribution of Ay and its uncertainty from 
z < 1 hosts (dark blue) and 1 < z < 2 hosts (olive). The red 
dashed line shows z > 2 hosts (error region not shown). 


= 1.31 mag, Anrf = 1.381 mag for Ay 
decrement of Ha over HjS follows 


Ay — 


2.5 ■ log(rQ]3s/rtheo) 


Au„ — A 


1 mag. For the Balmer 


(5) 




Fig. 10; [Oiii]/[Oii] as a function of (from top to bottom); red- 
shift, galaxy Ay, and SFR. The color of individual data points are 
according to the color bars. In all panels, solid lines denote the 
most-likely linear regression, with regions of 68 % probability in 
shaded areas. Dusty/normal GRBs are indicated by stars/circles. 


stellar mass M*, w here Ay =1 mag corresponds roughly to 
M* ~ lO*** Mq (e.g., Zahid et al.'2014 1 . 


The [O in]/[0 ii] ratio thus increases with redshift, whereas 
it decreases with Ay (or stellar mass), plausibly an effect of the 
efficient absorption of ionizing photons on dust. The net effect 
makes [O iii]/[0 ii] independent of SFR in our sample. 


5.3. Dust properties - distribution and redshift evoiution of 
gaiaxy Ay 

We measure the rest-frame attenuation of light from dust through 
comparing the observed Balmer decrement robs with the the¬ 
oretical expectation rtheo for a given electron temperature. To 
exploit the broad wavelength coverage and detection of more 
than two Balmer lines, we simultaneously minimize the various 
expected ratios against Ay for the four strongest Balmer lines. 



four lines correspond to Aho = 0.79 mag. Ah/s = 1.16 mag. Any 


'^We choose a MW-type extinction law to compare directly with 
other studies, and note that the difference to extinction laws from the 
SMC or LMC is small in the wavelength range of Hh to Ha. 


and similarly for the other ratios. In the calculation of Ay, we use 
a prior of Ay > 0 mag, and provide the median of the a-posteriori 
probability and Icr equivalent parameter range as error to repre¬ 
sent its distribution. In total, we constrain Ay for 71 events in the 
sample (Table [^, the distribution of which is shown in Fig. 11 
After applying the statistical analysis and cuts of 
Sects. |2.3.3| and \2.4\ the median 

\ med _ A q+0.4 , 


of our sample is 
Amed _ 0 . 8 ^q 3 mag and seems to be somewhat evolving 
with redshift (Fig. 111 . At redshifts lower than z < 1, there is 
a lower number of Ay > 1.5 mag hosts when compared to the 


z > 1 sample, and we measure A™‘^ = 0.55^[! (o mag. This value 


;+0.24 , 


is full y consistent with the value reported in Savaglio et al. 
( |2009 I. A™"* increases to A™‘^ = 1.2 + 0.3 mag at 1 < z < 2. It 
is A™'^ = 2 at z > 2 , but we note that the galaxies are 

significantly brighter than comparison samples of GRB hosts in 
this redshift range. 

The significance of this evolution in Ay with redshift how¬ 
ever is relatively modest; comparing the two distributions cen¬ 
tered around z ~ 0.6 (z < 1) and z ~ 1.5 (1 <z<2) with each 
other (Fig. the null-hypothesis that the two distributions are 
drawn from the same parent sample can be rejected with a sig¬ 
nificance of 2cr. At even higher redshifts further selection effects 
play an important role as well; high Ay hosts are not appearing 
because of the decreasing sensitivity SFR at high Ay (Fig.|^. 

The relative sparsity of high Ay events at low redshift is how¬ 
ever not easily explained by selection effects. At the average 
mass of the SF-weighted galaxy population of M* ~ 10^° Mq 
at z ~ 0.6 ([Perley ^ al.||2013b Hunt et ^[2014 1 , the average 


reddening for galaxies is Es-y ~ 0.35 mag (|Zahid et al.||2014 


and references therein). For GRBs to trace star formation rep¬ 
resentatively, half of the sample would be expected above this 
value and half below, while we observe (82 + 8 )% of all hosts 
below Eg-y = 0.35 mag at z < 1. 
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SFR(M0yi-i) 

Fig. 12: In the top panel, the blue solid histogram with 
error bars shows the SFR distribution and its uncertainty 
of z < 2 GRBs, where we have applied the correction 

for the over-proportionality of dusty GRBs in the sample 
(Sects. |2.3.3| and |2.4| l. Black data show all raw measurements 
without corrections. The lower panel shows cumulative distribu¬ 
tions of SFR and its uncertainty of z < 0.7 hosts (dark blue), 
0.7 < z < 1.3 hosts (olive) and, 1.3 < z < 2.0 hosts (red), in the 
same way corrected as the blue histogram. The black dashed line 
denotes hosts at z > 2. 


Although there are indications of an increasing dust-content 
in GRB hosts with redshift, the low significance prevents us 
from reaching stronger conclusions. Interestingly, however, a 
very similar trend - highest visual extinctions at z ~ 1.5 - 2 - 
exists in the complete sample of Ay values derived from after¬ 
glow data ( |Covino et al.||2013] l and in the evolution of the far 
UV attenuation in field galaxies ( Cucciati et al.|20r^ . It is also 
consistent with the somewhat higher fraction of dusty GRBs at 
higher redshift (Sect. [Z33| ). 


5.4. Star formation rates 

As introduced in Sect. |4.4| we can use the flux of nebular lines to 
derive constraints on the current star formation rate (SFR). Ha, 
in particular, is a common tracer of SF on short timescales of a 
few million years, and its dependence on the uncertainty in the 
dust-obscuration is only modest. 

All SFRs are given in Table and are corrected for the mea¬ 
sured Ay if available (Sect. or derived from scaling relations 
that intrinsically include the effect of dust (Sect. |4.4| l. Here we 
also add the systematic uncertainty of the flux calibration (Ta¬ 
bles and 1^. In the majority of all cases, however, the error 
budget is dominated by the uncertainty in the dust correction. 

Generally, all of the hosts in the sample are forming stars at 
a level of at least 0.02 Mgyr^*. Even though we are sensitive to 
galaxies with a very low SFR at z < 0.5, we detect none (but 
see |Rossi et aL]|2014[ for a quiescent galaxy at the position of 
GRB 050219A). Figure \T2\ shows the distribution of SFRs for 
89 hosts, out of which 68 are based on Ha, ten on H/5, eight on 
[On], and three on [Out]. The distribution is very broad, span¬ 
ning nearly four orders of magnitude between 0.05 Mq yr“^ up 
to around 200Moyr'' and peaks at 20Moyr“'. A large num¬ 
ber of galaxies with high SFRs either host dusty GRBs or are 
at z > 2 (where selection effects are important). The peak of the 


corrected SFR distribution for z < 2 GRBs (blue lines in Fig. 12 1 
is therefore an order of magnitude lower (2 Mq yr"^). 
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Fig. 13: SFR as a function of redshift. Blue circles and red stars 
represent individual galaxies. Large or small pentagons are the 
median values binned in redshift with bootstrapped errors at 
z < 2 and z > 2, respectively. Here, we have applied the cor¬ 
rection for the over-proportionality of dusty GRBs in the sam¬ 
ple (Sects. [T.3.3| and [Z^ . The gray shaded areas represent the 
average sensitivity as detailed in Sect. 3.4 using Ay = 1 mag 
(Sect. 1^. The three dashed lines show typical SFRs for differ¬ 
ently massive galaxies from the main sequence ( |Speagle et al.| 
2014[l for a qualitative estimate of the evolution of field galaxies. 


5.4.1. Redshift evolution of star formation rates 


As shown in Figs.[^and ^ there is clear evolution of the typi¬ 
cal SFR of GRB hosts with redshift. While at z ~ 0.5, the median 
SFR of a GRB host is ~ 0.6 Mq yr * (see also Savaglio et al. 


2009 1 , it rapidly increases by a factor of 25 to ~15Moyr ‘ at 


z ~ 2. Above z > 2, no further strong rise is evident in the data, 
and the median SFR levels out at around ~2QMq yr '. 

This behavior changes only slightly if we take into account 
that our selection might miss the faintest ~2Q% of hosts (total of 
~ 5) at 1 < z < 2 (Sect. 2.3 and Fig.|^, or that we have initially 
excluded four spectra of galaxies that could potentially be in this 
redshift range (Sect. in). 

The 30-40% quantile of SFRs at z > 1.3 is still an order of 
magnitude larger than the median SFR at low redshifts (Fig. [T^ 
lower panel). To be consistent at the 90% confidence level with 
no evolution in SFR, we would need to have missed by chance 
more than 40 GRBs (i.e., twice the sample size in the 1.3 < z < 
2.0 range) that are hosted by galaxies with SFR <2 Mq yr^'. 

At least at z > 2.5 there is an obvious effect of sensitiv¬ 
ity (shaded area in Fig. [T3|: only the most star-forming hosts 
(>20Moyr“^) have emission lines bright enough to be picked 
up by X-Shooter. We hence sample only the brightest end of 
the SFR distribution, and the comparison to the z < 2 data be¬ 
comes less meaningful. It is however clear that there is no further 
order-of-magnitude increase in the median SFR of GRB-selected 
galaxies above z > 2. 

The strong evolution in the SFR is not particularly surpris¬ 
ing: indeed it is similar to the redshift dependence of the UV lu¬ 
minosity of normal galaxies parametrized with e.g., the normal¬ 
ization of the main sequence (e.g., |Wuyts et al.|2011 [Whitaker 
|et al.|2012| l. To illustrate the evolution of field galaxies, we also 
plot redshift tracks for galaxies on the main sequence with stel¬ 
lar masses between M* = 10^ Mq and 10’' Mq in dashed lines 
in Fig.[T^ They indicate qualitative (but not quantitative) agree- 
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Fig. 14; SFR-weighted Ha-luminosity function at z = 0.8 
( Ly et aL]|2011| l in red compared to the histogram of SFRs 
from GRB hosts. The blue-solid histogram with error bars 
shows the SFR distribution and its uncertainty of our 0.55 < 
z < 1.05 GRBs, where we have applied the correction 

for the over-proportionality of dusty GRBs in the sample 
(Sects. |T.3.3| and [Z4| ). Black data are raw measurements in the 
respective redshift range without corrections. Shown curves are 
normalized to an area of 1. 



Fig. 15: In the top panel, the blue histogram with error bars 
shows the distribution of line broadenings (cr) and its uncertainty 
of z < 2 GRBs (Sects. [273.3| and [2^ . The black histogram shows 
raw measurements without corrections. In the bottom panel we 
plot cumulative distributions of cr, where black/blue/red colors 
denote redshift intervals centered around of z ~ 0.5, z ~ 1.1, 
and z ~ 1.6, respectively. To enhance clarity, the error regions 
around the z ~ 1.1 and z > 2 histogram are not shown. 


ment in the relative rise in SFR with redshift between the differ¬ 
ent samples. There are significant differences, in particular be¬ 
cause of a relative lack of galaxies in the high SFR regime at low 
redshift. 


5.4.2. Sparsity of high SFR GRB hosts at z < 1 

We illustrate the sparsity of high SFR galaxies hosting GRBs at 
low redshifts in Fig. [fusing the GRB host sample between z = 
0.55 and z = 1.05 ((z) = 0.81) and the dust-corrected luminosity 
function (LF) of Ha from Ly et al. ( 201 Ij ) at a similar redshift 
(Schechterparameters a = 1.6 and log(L*/(erg s“')) = 43.0). We 
multiply the LF by the SFR (or, in this case, Ha luminosity) since 
GRBs are expected to follow total SFR and not galaxy numbers. 

The SFR-weighted Ha-LF then represents the distribution 
that would be obtained from GRB hosts if they directly traced all 
star formation. Figure shows, however, a strong discrepancy 
between both distributions (see also |Boissier et al.|2013| l with the 
GRB hosts peaking at around an order of magnitude lower SFRs 
(~ IMoyr^*) than the SFR-weighted Ha-LF (~ lOMoyr^^). 
There is no observational bias in the sample that could explain 
this discrepancy. 


The virial or dynamical mass Md of a galaxy is then a func¬ 
tion of galaxy geometry, velocity dispersion era, and typical 
physical extent, or h alf-light radius r, of the galaxy via the virial 
theorem (see, e.g., Pettini et al. 2001 \ as Md = Ca-^^r^jG). 
C is related to the galaxy's geometry, density, and velocity 
anisotropy, with C = 5 in the case of an isotropic sp here (jPet- 
tini et al.|2001 1, or C ~ 3.4 for a disk-like geometry (Erb et al. 


2006bl l. 


No measurements of the effective radius for the bulk of the 
sample is available to us because they are spatially unresolved 
from the ground. In addition, a systematic uncertainty is present 
at high redshift because the emission line shape is dominated 
by luminous Hn regions, probing only a part of the gravita¬ 
tional potential. Despite all these uncertainties, the simple ve¬ 
locity broadening is an important tracer of galaxy properties as 
demonstrated in Figs. [T6| and [20] 


5.5. Line broadening 

The broadening of the nebular lines traces the turbulent and ro¬ 
tational motion of the hot gas within the gravitational potential 
of the host galaxy. The line-widths are well resolved in all but a 
few low-redshift cases because of the medium resolving power 
(4000 < R < 10000) of X-Shooter. The line-broadening from 
long-slit spectroscopy correlates reasonably well with velocity 
dispersions measured in spatially resolved data ( [Richard et al.j 
|2011| l, so we directly use the line-width cr as a proxy for the gas 
velocity dispersion crj. For the bulk of our targets, no evidence 
for rotation is seen in the spectra, but we caution that the ro¬ 
tational component contributes to the measured line-broadening. 


The peak of the line-width distribution (Fig. 151 at 40 km s“' 
suggests a dynamical mass on the order of Md ~ the 

lowest measured values (20kms“') imply Md ~ IO^'^Mq. At 
the high cr end, there are only three hosts above cr = 150 km s"\ 
namely GRB 070521 (z ~ 2 and with a large error on cr), 
GRB 090201 (z ~ 2), and GRB 120118B (z ~ 3) which indi¬ 
cate high dynamical masses (Md ~ lO" Mq). These three galax¬ 
ies are selected through dust-extinguished GRBs and their cr is 
comparable to those of the most luminous galaxies at the respec¬ 
tive redshift interval (e.g., [Erb et al.|200^ . 

With respect to redshift (Eig. 161, there are indications for 
an evolution such that lower redshift GRB hosts have - on aver¬ 
age - a smaller line-broadening. In particular, there is a lack of 
high cr galaxies at z < 1 (Eig. [T^ lower panel), while there is 
no obvious bias that would prevent us from selecting or observ¬ 
ing these galaxies in the first place. The average cr increases - 
qualitatively similar to the evolution of Ay in Sect. 5.3 - from 
cTmed = 40kms“' at z = 0.5 to 60kms“' < cr < 70kms“' at 
z > 1.0. A mild correlation {p - 0.02) exists also between Ay 
and cr, two quantities which relate to the total mass of a galaxy. 


in particular at higher stellar masses (e.g., Eorster Schreiber et al. 
[2O09t . 


More obvious in Eig. 16 is the strong correlation between 
SFR and cr (< 10“^ of all samples have no correlation). A simi- 
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Fig. 16: Line broadening cr as a function of (from top to bottom): 
redshift, Ay, and cr. In the top panel, large/small pentagons show 
the median values binned in redshift with bootstrapped errors 
at z < 2 and z > 2, respectively. The color coding of events 
corresponds to redshift, with the scale given by the color bar. In 
the lowest two panels, solid lines denote the most-likely linear 
regression, with error regions (68 % probability) in shaded areas. 


lar relation between luminosity and velocity dispersion (L oc cr^) 
is well-established in the local Universe. It was first observed in 
ellipticals by Faber & Jackson ( 1976[), a nd later also in z ~ 2 
star-forming galaxies ( |Erb et ah 2006b| l. Murray et al. (2005 i 
describe this correlation as a consequence of moderation of star 
formation activity through momentum-driven winds from super¬ 
novae together with radiation pressure. Star formation would 
drive-out gas above a critical luminosity Lm AfgCjG ■ cr^, 
where fg is the gas fraction. At the bolometric luminosity Lm 
a star-forming galaxy then self-regulates its SFR and reaches its 
maximum luminosity. 

The correlation between SFR and cr in Fig. is consistent 
with S FR oc cr^, but there is considerable dispersion in the data. 
To test for a possible redshift evolution, we fit galaxies below 
and above z = 1 with a fixed proportionality of cr"'. This returns 
within errors similar normalizations with minor evidence (1.5 cr) 
for the higher-redshift galaxies shifted towards higher SFRs for 
a given cr. This would indicate higher gas-fractions in the model 
of |Murray et al. ( 2005| l. 


5.6. Metallicity 

5.6.1. Metallicity measurements from strong-line diagnostics 

One of the most fundamental - but also most difficult to mea¬ 
sure - characteristic of galaxies is their metal abundance. While 
for GRB afterglows, metallicity can be probed directly and ac¬ 
curately in absorption at z > 1.8 through their DLAs (e.g.. 
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Vreeswijk et al. |2004[ |Savaglio| 2006[ [Prochaska et al.|[2007t 


avoidable to resort to diagnostic ratios of strong emission lines 

(e.g.,|Erb et al.|2006a 

[Maiolino et al.[[2008[|. These strong line 

diagnostics use metal 

icity-dependent ratios of e.g., [On] and 

[Oiii] to H/J (e.g., Pagel et al.[[19791 McGaugh[[ 19911 Zaritsky 

et al.|[1994|l, [Nn] to Ha ([Alloin et al. 

|1979[ 

Pettini & Pagel 


2004| l, or [Nii] to [On] (e.g., Kewley & Dopita|2002[ ) to derive 


an oxygen abundance 12-i-log(0/H). For an extensive discussion 
on the strong line ratios, we refer to |Kewley & Dopita| ( |2002| l. 

In particular at high redshift, a reliable determination of the 
metal content in the hot gas phase becomes challenging through 
these line ratios (e.g., |Maiolino et al.||20081 l. As we have seen 
in the previous paragraphs, the physical conditions, e.g., ion¬ 
ization parameter, SFR, possibly dust content, and gas-fractions 
are different from those of low-redshift galaxies on which most 
of the various strong-line diagnostics are calibrated. Currently, 
only very few directly measured abundances from high-redshift 
galaxies are available (e.g., [Christensen et al.|2012) l and none of 
the strong-line ratios are thus validated at z ~ 2. Any metallicity 
that is derived from these strong-line ratios must thus be taken 
cautiously (e.g., |Steidel et al.|2014| ). 

In addition, we are faced here with the problem that our sam¬ 
ple spans a broad range of redshifts, broad range of SFRs and 
host masses, and a different set of lines is detected for every 
host. Because of the different calibrations inherent in different 
strong-line diagnostics, any comparison of them needs thorough 
consideration (e.g., Kewley & Ellison|2008| l. 

Given the aforementioned constraints, we decide to base our 
metallicity measurements on the single set of calibrators from 
Nagao et al. (2006i and Maiolino et al. ( |2008| l, which provide 


a homogeneous set of line diagnostics for different line ratios. 
Specifically, we use the five diagnostics based on R 23 , [On] and 
[Neiii], [Out] and [Nii], Ha and [Nii], [Nii]and [On] and si¬ 
multaneously minimize 12-i-log(0/H) against the available data. 
Via the line-flux measurements and errors (i.e., likelihood dis¬ 
tributions of fluxes), we construct probability density functions 
(PDFs) for each diagnostic ratio, as well as one PDF for the com¬ 
bination of them. Sectionj^and Fig. jB.lj in the Appendix provide 
more details on the procedure and individual events. 

We also provide all line-flux measurements and errors in Ta- 
bles|^and[^ so oxygen abundances can be reproduced and recal¬ 
culated on any given metallicity scale with the preferred strong¬ 
line diagnostic. 


5.6.2. Sample properties 

In total, we measure well-constrained gas-phase metallicities for 
44 events (Table in a redshift range between z = 0.3 and 
z = 3.4 with a median of Zmed = LO. The oxygen abundances 
are distributed between 12 + log(0/H) = 7.9 (GRB 120118B) 
and 12 -h log(0/H) = 9.0 (GRB 051117B) with a peak of the 
distribution at 12 - 1 - log(0/H) ~ 8.5 (Figure [TT]). 

Intriguingly, we detect a non-negligible fraction of GRB 
hosts at high metallicities (Z > Zq, Fig. see also [Levesqu^ 
jet al.|2010ct [Elliott et al.||2013| l. It is thus obvious that there is 
no strict cut-off in host metallicities at least up to the solar value. 
We thus directly confirm the predictions from previous studies 
which indicated that GRBs form in various types of star-forming 
galaxies, even metal rich ones ([Berger et al.|2007[ [Kriihler et al.[ 
[20TT]|Perleyetal.|2013b|[2015c|l. 

The number of hosts with a metallicity above solar, however, 
is relatively small. From the sub-sample of galaxies with metal¬ 
licity measurements, we derive fractions of 24+ 10%, 26+ 13%, 
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Fig. 17; The top panel shows the histogram of measured oxygen 
abundances (12 + log(0/H)) derived through strong-line diag¬ 
nostic ratios in blue following Sects. [Z331 and [Z4] and in black 
for all data. In the bottom panel we plot cumulative distribu¬ 
tions of the GRB host metallicities in different redshift intervals, 
where blue/red/olive colors denote redshift intervals below z - I 
{Zmed = 0.6), 1 < z < 2 (Zmed = 1-4), and 2 < z (Zmed = 2.4), 
respectively. Zq is indicated by the vertical, gray dashed line. 
Error-regions of the 1 < z < 2 sample have a similar size to the 
one for z < 1, but are not shown to enhance clarity. 


7^]“% for z< 1,1 <z<2, and 2 < z. Under the assump¬ 
tion that the emission-line width cr traces oxygen abundance 
as shown in Sect. |5.6.5[ we can estimate metallicities also for 
those galaxies for which the strong-line diagnostics can not be 
applied. For instance, the host of GRB 050525A at z = 0.6 


is faint {R - 25.7 mag. 

Hjorth et al. 

2012|l, has a narrow line 

width and low SFR, anc 

is thus like 

y metal-poor (the corre- 


including all galaxies of this work, the fractions of super-solar 
metallicity hosts become 16 + 7%, 25 + 9%, 4^2*^ at z < 1, 
1 < z < 2, and 2 < z, respectively. 

At the highest redshifts z > 3, emission-line metallicities are 
traditionally very expensive and challenging to obtain. We derive 
through observations of [Neiii]for GRBs 110818A (z = 3.36), 
111123A (z = 3.15), and 120118B (z = 2.94) metallicities of 
12 H- log(0/H) = 8.25+“'^’, 12 -H log(0/H) = 8.01 + 0.28, and 
12 -H log(0/H) = 7.89 ^q? 2, respectively. Comparing this with 


results on galaxy metallicities at z ~ 3.5 and above (Maiolino 
let al.|2008t|Laskar et al.|20rT] l, we see that those three events are 
probing the highest measured metallicities of LBGs at the given 
redshift. This is consistent with the significant line broadening 
and high SFRs (Sects. 5.5 and |5.4| ). 


5.6.3. Evolution with redshift 

To test for an evolution of the metal content in the gas-phase 
with redshift, we divide the galaxies into low (z < 1, Zmed = 0.6), 
medium (1 < z < 2, Zmed = 1-4), and high redshift (z > 2, 
Zmed = 2.4) samples and plot histograms of 12 + log(0/H) in the 
bottom panel of Fig. Remarkably, the cumulative distribu¬ 
tions appear similar in the two low-redshift intervals. In contrast, 
the z > 2 sample is somewhat shifted towards lower metallici¬ 
ties. 



Fig. 18; The blue-solid histogram with error bars shows 
the metallicity distribution and its uncertainty of 29 0.3 < 
z < 1 GRB hosts ((z) = 0.65). We applied a correc¬ 

tion for the over-proportionality of dusty GRBs in the sample 
(Sects. 2.3.3| and |2.4| ) and for completeness from the line broad¬ 
ening. Black data show raw measurements. The red line is the 
expectation from the stellar-mass function at a similar redshift 
( |Tomczak et al.|[2014| l if GRBs traced star formation without 
metallicity dependence. Shown curves are normalized to an area 
of 1. 


Cosmic-chemical evolution arguably proceeds most rapidly 
at higher redshifts. Between z ~ 0 and z ~ 0.7 the decrease 
in 12 H- log(0/H) for a galaxy of a given mass seems relatively 
modest at around 0.1-0.2 dex ( |Tremonti et al. |2004[ |Savaglio| 
et al.|2003 1, even though this difference in redshift corresponds 


to over 5 Gyr in the evolution of the Universe. The displacement 
of the mass-metallicity relation at high redshift is likely much 
larger than this, reaching ~0.5 dex at the high mass (10/* Mn) 


end, and 1 dex at 10® Mq at z ~ 3.5 (Maiolino et al. 2008 Tron 


|coso et al.||2014| l. This leads to relatively low gas-phase metal¬ 
licities in the sample although we are only sensitive to the more 
luminous GRB hosts at z > 2. 


5.6.4. Metallicity dependence of GRB hosts 


To understand the effects of metallicity on GRB progenitors 
and host galaxies and the implications for their ability to trace 
star formation, we summarize the primary results from the pre¬ 
vious sections; First, GRB hosts with oxygen abundances of 
12 + log(0/H) ~ 8.5 are common in our sample. Second, the 
z < 1 and 1 < z < 2 GRB host metallicity distributions appear to 
be very similar. Third, there is a modest fraction (~ 20% at z < 1) 
of GRBs in Z > Zq galaxies. And last, the maximum metallicity 
of the most luminous GRB hosts at z ~ 3 is 12 -i- log(0/H) ~ 8.4. 


In particular the relatively low fraction of GRB hosts with 
very high oxygen abundances points strongly towards a metal¬ 
licity dependence in low-z GRBs. In Fig. 


18 


we compare the 
metallicity distribution derived from 29 0.3 < z < 1 GRB hosts 
((z) = 0.6) with the expectation from the contribution of field 
galaxies to the global SFR density at a similar redshift. 

As no sufficiently complete sample of field-galaxy metallic¬ 
ities is available in the literature at z ~ 0.6, this comparison is 
based on the stellar-mass function from deep photometric sur¬ 


veys. To derive the red lines in Fig. 18 we start with the double 
Schechter parameterization of the galaxy-mass function for Star- 
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forming galaxies at 0.5 < z < 0.75 from Tomczak et al. (2014 
Passively evolving galaxies are hence implicitly excluded from 
this comparison. Using the mean-value of the SFR-stellar mass 
relation and its dispersiorf^ we then calculate the 50% quartile 
of SFR for a given stellar mass. 

Finally, we use the parametrization of |Mannucci et al.| ( |2010) l 
to convert the SFR, stellar mass pairs into a metallicity in the 
same scale as our GRB measurements and weight the resulting 
galaxy-metallicity function with the SFRp] 

The fraction of star formation in galaxies above Zq can then 
be obtained via the integr al under the red curve of Fig. and is 
around 50 Similarly, Campisi et al. (20111 use N-body sim¬ 
ulations with semi-analytical galaxy-formation models to predict 
around ~ 50% of GRBs in galaxies of super-solar metallicity in 
cases where GRBs traced all star formation equally (their Fig. 

4). 

Our results therefore directly confirm studies that at low red- 
shift, GRB host properties are affected by a tendency of GRBs 
to occur in environments with a metallicity below the solar value 
(e.g., Floc’h et al.|2006[|Kocevski et al.|2009)[Levesque et ah 


2010b| Graham &; Fruchter|2013 |Perley et al.||20 

3b |Boissier 

et al.|2(J13| Vergani et al.|2015l Schulze et al. 2015 

Perley et al. 


2015c I. The scarcity of high Ay and high SFR galaxies as dis¬ 
cussed in previous sections further strengthens this interpreta¬ 
tion. 


The physical driver behind the environmental bias in the 
GRB production efficiency could in principle relate to either (or 
both), metallicity or a high stellar mass. Because of the well- 
studied impact that metallicity has on stellar evolution, and the 
lack of an obvious link between GRB formation and a galaxy’s 
stellar mass, we consider a metallicity effect as the most logical 
explanation for our observations. This interpretation is further 
supported by the increasing stellar mass of GRB hosts with in¬ 
creasing redshift ( |Perley et al.|2015c] l. 

The observed metallicity dependence is unlikely to play a 
dominant role at z > 3. Already z > 2, the metallicity distri¬ 
bution is somewhat shifted towards lower metallicities although 
the selection biases against low-luminosity hosts become more 
severe. At z ~ 3, galaxies of solar metallicity are scarce even 
among the most luminous galaxies ( jTroncoso et al.|[2014| l. Be¬ 
cause the metallicity of star-forming galaxies in general is de¬ 
creasing, GRBs are expected to form also more frequently in 
the most luminous galaxies. Therefore, we expect GRB hosts at 
high redshift to represent star-forming galaxies much better than 
at z < 1 (see also|Greiner et al.|2015a[|Perley et al.|2015c[). 




Fig. 19; Gas-phase metallicity as a function of (from top to bot¬ 
tom); redshift. Ay, SFR, and cr. Black data and error bars rep¬ 
resent individual galaxies. The color coding corresponds to red¬ 
shift, with the scale given by the color bar. Zq is indicated by the 
horizontal, gray dashed line. 


There is already evidence that the implications of the metal¬ 
licity dependence are weaker at z ~ 2 ( |Perley et al.||201 Jbj 
Schulze et al.|2015] l, than in z < 1 GRB-selected galaxies. 


''Only subtle differences are introduced in the result when using the 


mass 

ilar n 

15 

2014 

-function of star-forming galaxies from|Ilbert et al.H2 
idshift interval. 

dere, we assume an intrinsic dispersion of 0.3 dex ( 

1 and a slope of 0.8 in the stellar-mass range b 

013| in a sim- 

Speagle et al. 
slow IO'^Mq 

(|Dunne et al. 

120091 

Santini et al.||2009| |Sobral et al.||20141( and 0.55 

above (IWhita 

cer et a 

1. 20121. 


the host’s 12 log(0/H) and SFR to derive a stellar mass via|Mannucci| 


et al. 

tr 


^j2010| to compare it with the SFR-weighted stellar-mass function. 
Similar values are inferred through the analysis of|Perley et al. 


( |2013b^ , jFlunt et ^(|2014| and|Sobral et al!](|2014^, which use deep 
photometric surveys ( [Kajisawa et al.|2009[|Ilbert et al.|2013[ ( to calcu¬ 
late a galaxy stellar mass of Mq at z ~ 0.6 above or below 

which half of all star formation occurs. This implies a metallicity be¬ 
tween Z = 0.9 Zq and Z = 1.4Zq in our oxygen-abundance scale jM an- 1 
nucci et al. 2010|. 


5.6.5. Correlations of metallicity with other host properties 

Figure shows the dependence of host metallicity with red¬ 
shift, dust reddening, SFR, and emission-line width. A correla¬ 
tion exist (p = 0.0007 of all iterations return no dependence) 
between Ay and oxygen abundance similar to 12 -i- log(0/H) = 
8.36 + 0.17 ■ Ay/(mag). This relation is driven by z < 2 events, 
as only a few z > 2 hosts have sufficient data to constrain both 
Ay and 12-i-log(0/H) at the same time. At somewhat higher red¬ 
shift, a correlation between metallicity and dust content is also 
seen directly along GRB sight-lines as traced by the GRB after¬ 
glow and its DLA ( |Zafar & Watson|2013| l. 

There are no strong correlations between SFR or cr and 
12-1- log(0/H) in the full host sample, but this is the result of the 
large redshift range present in Fig. 19 In Fig. 20 we plot oxy¬ 
gen abundance versus Ay, SFR, and cr divided into three red¬ 
shift bins, Z < 1 with Zmed ~ 0.6, 1 < Z < 2 with Zmed ~ 1 -4, 
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Fig. 20; Gas-phase metallicity as a function of z and (from top to 
bottom): Ay, SFR, and cr as well as z < 1, 1 < z < 2, and 2 < z 
from left to right. Solid lines denote the best linear regression, 
with regions of 68 % probability in shaded areas. The last row 
shows also the relations from the lower-redshift panels in dashed 
lines. Zq is indicated by the horizontal, gray dashed line. 


and 2 < z with Zmed ~ 2.4. A mild correlation exists between 
SFR and 12 -I- log(0/H) similar to 12 H- log(0/H) oc 5 iz/; 0 - 20 ± 0.06 
(p - 0.004). Using a common slope, but a different normaliza¬ 
tion in the three redshift ranges, we fit all data simultaneously 
and derive 


12 -H log(0/H)2<i = 8.63 H- 0.20 ■ logC^FR/lOMoyr"') 

12 -H log( 0 /H)i <^<2 = 8.49 H- 0.20 ■ log(5FR/10Mo yr"') 

12 -H log( 0 /H) 2 <j = 8.20 H- 0.20 ■ log(5FR/10Mo yr"') (g) 

with uncertainties of 0.05 in the intercept in the earlier two, 
and 0.11 in the z > 2 case. The root mean square (RMS) dis¬ 
persion in the z < 1 and 1 < z < 2 redshift range is 0.25 dex 
each. 

The decrease in log(0/H) for a given SFR between z ~ 0.6 
and z ~ 1.4 is 0.14 + 0.08 dex, but 0.4 + 0.1 dex for the z > 
2 sample (Fig. [20l i. This is consistent with the relatively slow 
evolution of the mass-metallicity relation at lower redshift, as 
well as a more rapid change at z > 1.5. 

The tightest and most signihcant correlation {p < 0.0001) 
is present between 12 + log(0/H) and cr as 12 + log(0/H) oc 
^o.64±o.i6 at z < 2. We again ht all data together with a common 


slope and free normalization in the three redshift intervals and 
obtain for cr 

12 + log(0/H),<i = 8.70 + 0.64 ■ log(o-/100kms-') 

12 + log( 0 /H)i<j <2 = 8.66 + 0.64 ■ log(o-/100kms^‘) 

12 + log( 0 /H) 2 <i: = 8.31 + 0.64 ■ log(o-/100kms-') (7) 

with uncertainties of 0.05/0.06 in the intercept in the earlier 
two, and 0.12 in the z > 2 case. The behavior is qualitatively 
very similar to relation between SFR and oxygen abundance but 
with a lower scatter, i.e., a tighter correlation. 

Low- and medium-redshift samples return comparable best- 
ht relations. In both cases cr = 100km s“^ roughly corresponds 
to a metallicity of Z a; Zq. At z > 2, this value shifts downward 
by 0.4 + 0.1 dex to an oxygen abundance of 12 + log(0/H) a; 8.3, 
but is affected by large measurement uncertainties. The RMS 
dispersion in the z < 1 and 1 < z < 2 redshift range is only 0.19 
dex and 0.20 dex, respectively. Given that the typical error in the 
metallicity measurement via the strong-line diagnostic is already 
0.11/0.16 dex, this implies that the very simple measurement of 
cr provides already a good estimate of the gas-phase metallicity. 

Taken at face value, there is no statistically significant cor¬ 
relation in the highest-redshift bin (z > 2) between SFR and 
12 + log(0/H) as well as cr and 12 + log(0/FI). Because of the 
large measurement uncertainties however, the data are consistent 
within errors with the lower-redshift slope, but, as argued above, 
a lower normalization. 


6. Summary 

We present data and initial results from X-Shooter emission¬ 
line spectroscopy of star-forming galaxies selected through long 
GRBs in a redshift range 0.1 < z < 3.6. In total, we discuss 
observations of 96 individual targets which allow us to trace 
the physical properties of GRB hosts through cosmic time. The 
galaxies have a median redshift of Zmed = 1.6 and most of them 
are in the range 0.5 < z < 1.5 (35%) and 1.5 < z < 2.5 (41%). 
We carefully test for selection effects in the sample, and statis¬ 
tically correct for a somewhat higher fraction of dusty GRBs 
(35 + 5%) when compared to the one of representative GRB 
samples (20 - 30 %). 

Our observation yield 16 new GRB redshifts, which increase 
the redshift completeness of previous samples to 87% (TOUGH) 
and 97% (BAT 6 ), with median redshifts of Zmed = 2.06 + 0.18 
and Zmed = 1.67 + 0.15, respectively. 

After excluding a very nearby event (GRB 100316D), a pos¬ 
sibly short GRB (GRB 100816A), and five galaxies without clear 
emission lines we arrive at a total number of 89 emission-line 
spectra from galaxies selected by GRBs. Based on these data, we 
provide new scaling relations to derive SFRs from emission lines 
of [O It], Hj3, and [O iii] even in the absence of a reliable dust cor¬ 
rection (but with a large scatter). We measure galaxy-integrated 
star formation rates (SFRs), visual dust-attenuations (Ay), line 
broadening (cr), and oxygen abundances (12 + log(0/H)) for our 
targets (Table|^. This yields the largest and most comprehensive 
sample of emission-line spectroscopy of GRB hosts available to 
date and allows us to study the distribution of the physical pa¬ 
rameters and their evolution with redshift. 

Our main results can be summarized as follows: 

1. The distribution of Ay for long GRB hosts is broad with the 
bulk of the targets at Ay < 1.5 mag and median A™"* ~ 
0.75 mag at z < 2. We observe a mild trend with redshift 
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such that GRB hosts at z = 1.5 tend to have larger dust con¬ 
tent (A™"* ~ 1.2 mag) and line broadening (cTmed ~ 70 kms) 
on average than low-redshift hosts (A™*^ ~ 0.6 mag, and 
coined ~ 40 kms). The significance of each of these results is 
only around 2cr, but the evolution in Ay is remarkably con¬ 
sistent with a similar behavior observed for GRB afterglows 
( |Covino et al.|2013| l. 

2. There is a strong evolution of the average GRB host proper¬ 
ties with redshift. The median SFR evolves from S f/Jmed = 
0.6 Mq yr“’ at z ~ 0.6 up to S T’Rmed = 15 Mq yr“' at z ~ 2 
above which it does not increase significantly any further. 

3. The ratio of [Oiii]/[Oii] of the galaxies increases with red- 
shift as well, which moves the star-forming galaxies away 
from the location of local galaxies in the BPT diagram. The 
value changes from [Oiii]/[Oii] ~ 0.9 at z ~ 0.3 by a factor 
of three to [Oiii]/[Oii] ~ 3 at z ~ 3. While roughly similar 
at z ~ 0.6 and z ~ 1.4, the metallicity distribution of GRB 
hosts starts to shift to lower metallicities above z - 2. 

4. We find a strong correlation between cr and SFR broadly 
consistent with S FR oc cr^ as known from the local Faber- 
Jackson relation. There is no statistically significant evidence 
in the data that this proportionality evolves with redshift. 

5. Further correlations exist between oxygen abundance and vi¬ 
sual attenuation (12 + log(0/H) cx 0.17 ■ Ay), SFR (12 -H 
log(0/H) oc SFR^ ^), and a (12 H- log(0/H) oc In the 
last two cases, we observe a relatively similar behavior at 
z < 1 and 1 < z < 2, but an evolution towards lower metal¬ 
licities at z > 2 by ~ 0.4 dex for a given SFR or cr. 

6. We detect several hosts with high oxygen abundances which 
rule out a strict metallicity cutoff for GRB hosts below the 
solar values (in the adopted scale). At the same time, how¬ 
ever, the fraction of GRBs hosted in Z > Zq galaxies at 
z < 1 is smaller (18 + 7 %) than what would be expected 
from the contribution of similarly metal-rich galaxies to the 
total cosmic SFR (> 50%). A mechanisms is thus at place 
that quenches GRB formation at the highest metallicities. 

At z = 3, most of the star formation takes place in galaxies 
of Z < 0.5 ■ Zq. At the same time GRB hosts with similar metal¬ 
licities are relatively common in our sample at lower redshift. 
Under the conditions that the probability of producing GRBs be¬ 
low a certain metallicity is constant (and not a strong function 
of metallicity itself) then by z ~ 3 GRBs host will probe a large 
fraction of the total star formation. In absence of further sec¬ 
ondary environmental factors, GRB hosts would then provide 
an extensive picture of high-redshift, star-forming galaxies. This 
conclusion is similar to the one reached through studying metal- 


licities of GRB-DLAs (jFynbo et aL||2008 |Arabsalmani et al. 

2015^, the UV-Iuminosity function at a similar redshift Pakobs- 

son et al.|20051 ISchuIze et al.|2015 

Greiner et al.|2015all, cos- 

mological simulations (Chisari et al. 

20101 or modeling the mass 


distribution of GRB hosts ( |Kocevski et al.|2009| l. 


Our results can be understood in the cosmological context 


and a picture drawn by recent studies ( 

Modjaz et al.|2008l Ko- 

cevski et al.|2009 Boissier et al.|2013 

Perley et al.|2013b Ver- 

gani et al. 2015||Wang & Dai|2014[ Trenti et al.||2015[ Perley 

et al.|2015c Schulze et al.|2015|l in which the GRB host prop- 


erties at lower redshift (z < 1-2) are driven by the GRB’s ten¬ 
dency to occur in lower-metallicity galaxies without fully avoid¬ 
ing metal-rich ones. The scarcity of high metallicity, high Ay, 
and high SFR galaxies at low redshift strongly supports this in¬ 
terpretation. 
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Table 1: X-shooter spectroscopic observations 


GRB host 

Gal. Eg-Y 
(mag) 

Exposure time (s) 

UVB VIS NIR 

UVB 

Slit width 

VIS NIR 

Obs. date 

References 

GRB 050416A 

0.026 

4x900 

4x900 

12 X 300 

I'.'O 

0'.'9 

0r.'9 

2011-Jan-19 

(1), (2) 

GRB 050525A 

0.082 

4 X 1050 

4 X 1050 

12 X 350 

T.'O 

0'.'9 

(y:9JH 

2012-Sep-18 

(1) 



4x630 

4x664 

4x695 

T.'O 

0'.'9 

a:9JH 

2012-Sep-21 


GRB 050714B 

0.048 

12 X 940 

12 X 900 

36 X 325 

0'.'8 

0'.'7 

0'.'6 

2013-Mar-19 

(1) 

GRB 050819 

0.100 

2x900 

2x900 

4x450 

I'.'O 

0'.'9 

(y.'9 

2010-Oct-29 

(1), (3) 

GRB 050824 

0.030 

4x900 

4x900 

8 x450 

I'.'O 

0'.'9 

(y:9JH 

2011-Sep-06 

(1) 

GRB 050915A 

0.022 

6 X 1200 

6 X 1200 

12 X 600 

I'.'O 

0'.'9 

(y:9JH 

2011-Sep-06 

(1), (3), (4), (5) 

GRB 051001 

0.013 

2 X 1800 

2 X 1800 

4x900 

T.'O 

0'.'9 

(y:9JH 

2010-0ct-30 

(1), (3) 

GRB 051016B 

0.035 

2x600 

2x600 

4x300 

T.'O 

0'.'9 

(Y.'9 

2009-Dec-23 

(1) 

GRB 051022A 

0.052 

2 X 1200 

2 X 1200 

4x600 

T.'6 

T.'5 

a:9JH 

2012-NOV-13 

(6) 

GRB 051117B 

0.048 

1 x720 

1 x720 

2x360 

T.'O 

0'.'9 

or.'9 

2009-Dec-23 

(1) 

GRB 060204B 

0.015 

4x900 

4x900 

9x300 

T.'O 

0'.'9 

(Y.'9 

2015-Mar-05 


GRB 060306 

0.030 

4x 1800 

4 X 1800 

12 X 600 

T.'3 

T.'2 

T.'2 

2011-Sep-01 

(1), (3), (5) 

GRB 060604 

0.037 

4x 1800 

4 X 1800 

8x900 

T.'O 

0'.'9 

(Y.'9 

2010-0ct-30 

(1), (3) 

GRB 060707 

0.019 

4x860 

4x850 

12 X 300 

T.'3 

T.'2 

T.'2 

2009-Oct-22 

(1), (7) 

GRB 060719 

0.007 

4 X 1500 

4 X 1500 

12 X 500 

T.'O 

0'.'9 

or.'9 

2011-0ct-20 

(1), (3), (5) 

GRB 060729 

0.048 

4x810 

4x844 

4x875 

T.'O 

0'.'9 

(y:9JH 

2012-Sep-21 

(1) 



4x520 

4x554 

4x585 

T.'O 

0'.'9 

(y:9JH 

2012-Sep-21 


GRB 060805A 

0.021 

10 X 940 

10 X 900 

30 X 320 

0'.'8 

an 

0'.'6 

2013-Mar-19 

(1), (3) 

GRB 060814 

0.034 

2 X 1800 

2 X 1800 

6x600 

T.'O 

(Y.'9 

(Y.'9 

2010-Feb-17 

(1), (3), (5) 

GRB 060912A 

0.045 

2x900 

2x934 

2x965 

T.'3 

T.'2 

T.'2 

2012-Sep-21 

(1) 

GRB 060923B 

0.129 

2 X 1800 

2 X 1800 

6x600 

T.'O 

0'.'9 

a.'9JH 

2013-Mar-20 

(1) 

GRB 060926 

0.137 

4x 1800 

4 X 1800 

12 X 600 

T.'O 

0'.'9 

0r.'9 

2011-Apr-24 

(18) 

GRB 061021 

0.049 

12 X 940 

12 X 900 

36 X 320 

T.'O 

0'.'9 

(y:9JH 

2013-Mar-21 

(1) 

GRB 061110A 

0.076 

4x 1800 

4 X 1800 

12 X 600 

T.'O 

0'.'9 

(y:9JH 

2011-0ct-20 

(1) 

GRB 061202 

0.131 

4x900 

4x900 

9x300 

T.'O 

0'.'9 

(y.'9 

2014-NOV-29 


GRB 070103 

0.058 

2x900 

2x900 

4x450 

T.'O 

0'.'9 

(Y.'9 

2010-Oct-29 

(1), (3) 

GRB 070110 

0.013 

4 X 1680 

4 X 1680 

8x840 

T.'O 

0'.'9 

(Y.'9 

2010-Oct-29 

(1), (7) 

GRB 070129 

0.120 

4x900 

4x900 

12 X 300 

T.'O 

0'.'9 

(y.'9 

2009-Dec-23 

(1) 

GRB 070224 

0.049 

8 x940 

8 x900 

24 X 300 

T.'O 

0'.'9 

(Y.'9 

2013-Mar-20 

(1), (3) 



4x900 

4x920 

12 X 320 

T.'O 

0'.'9 

(Y.'9 

2013-Mar-20 


GRB 070306 

0.024 

6x900 

6x900 

6x900 

T.'O 

0'.'9 

(Y.'9 

2010-Dec-09 

(1), (5), (8), (9) 

GRB 070318 

0.015 

4x910 

4x885 

12 X 300 

T.'O 

0'.'9 

(Y.'9 

2010-0ct-06 

(1) 

GRB 070328 

0.031 

7x844 

7x810 

8x875 

T.'O 

0'.'9 

a:9JH 

2012-Sep-21 

(1) 



1 x589 

1 x585 


T.'O 

0'.'9 

(y:9JH 

2012-Sep-21 


GRB 070419B 

0.076 

4x 1800 

4 X 1800 

12 X 600 

T.'O 

0'.'9 

(Y.'9 

2011-0ct-20 

(1) 

GRB 070521 

0.023 

4x900 

4x900 

12 X 300 

T.'O 

0'.'9 

(Y.'9 

2013-Apr-03 

(5), (10) 



4x900 

4x900 

12 X 300 

T.'O 

0'.'9 

(Y.'9 

2013-May-07 


GRB 070802 

0.023 

4x 1800 

4 X 1800 

8x900 

T.'O 

0'.'9 

0r.'9 

2010-0ct-30 (1), (5), (9), (11) 

GRB 071021 

0.057 

4x 1800 

4 X 1800 

12 X 600 

T.'O 

0'.'9 

0r.'9 

2011-Sep-06 

(3), (5) 

GRB 080207 

0.020 

4x 1800 

4 X 1800 

12 X 600 

T.'O 

0'.'9 

(y.'9 

2011-Apr-24 

(3), (5), (10) 

GRB 080413B 

0.032 

3 X 1800 

3 X 1800 

6x900 

T.'O 

0'.'9 

(Y.'9 

2010-Aug-15 

(23) 

GRB 080602 

0.025 

4x900 

4x900 

12 X 300 

T.'O 

0'.'9 

(Y.'9 

2013-Aug-18 

(10) 



4x900 

4x900 

12 X 300 

T.'O 

0'.'9 

(Y.'9 

2013-Aug-19 


GRB 080605 

0.118 

4x885 

4x910 

12 X 300 

T.'O 

0'.'9 

0r.'9 

2011-Apr-26 

(9), (12) 

GRB 080804 

0.014 

4x 1350 

4 X 1350 

12 X 450 

T.'O 

0'.'9 

0r.'9 

2012-Apr-19 


GRB 080805 

0.037 

4x885 

4x910 

12 X 300 

T.'O 

0'.'9 

(y.'9 

2011-May-13 

(9) 

GRB 081109 

0.017 

2 X 1200 

2 X 1200 

4x600 

T.'6 

T.'5 

(y:9JH 

2012-Oct-13 

(5), (9) 



2 X 1200 

2 X 1200 

4x600 

T.'6 

T.'5 

a:9JH 

2012-NOV-14 


GRB 081210 

0.066 

4x900 

4x900 

12 X 300 

T.'O 

T.'O 

(y.'9 

2014-Oct-19 


GRB 081221 

0.019 

4x 1800 

4 X 1800 

12 X 600 

T.'3 

T.'2 

T.'2 

2011-Sep-01 

(5), (14) 

GRB 090113 

0.072 

2 X 1800 

2 X 1800 

4x900 

T.'O 

0'.'9 

a:9JH 

2011-Sep-06 

(3) 

GRB 090201 

0.059 

2 X 1200 

2 X 1200 

8 x300 

T.'O 

0'.'9 

or.'9 

2012-Dec-25 




2 X 1200 

2 X 1200 

8 x300 

T.'O 

0'.'9 

(y.'9 

2012-Dec-31 


GRB 090323 

0.021 

3 X 1560 

3 X 1560 

9x520 

T.'O 

0'.'9 

(y.'9 

2011-Mar-ll 

(15) 

GRB 090407 

0.059 

4x 1800 

4 X 1800 

8x900 

0'.'8 

0'.'7 

0'.'6 

2010-0ct-30 

(3), (5) 

GRB 090926B 

0.020 

2 X 1200 

2 X 1200 

8x300 

T.'O 

0'.'9 

a:9JH 

2012-NOV-14 

(9) 



2 X 1200 

2 X 1200 

8x300 

T.'O 

0'.'9 

a:9JH 

2012-Dec-17 


GRB 091018 

0.025 

4x600 

4x600 

4x600 

T.'O 

0'.'9 

(y.'9 

2009-0ct-19 

(16) 
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Table 1: X-shooter spectroscopic observations (continued) 


GRB host 

Gal. Eg-Y 
(mag) 

Exposure time (s) 

UVB VIS NIR 

UVB 

Slit width 

VIS NIR 

Obs. date 

References 

GRB 091127 

0.033 

4 X 1500 

4 X 1500 

8 x750 

I'.'O 

0'.'9 

0r.'9 

2009-Dec-02 

(17) 

GRB 

100316D 

0.101 

4 X 1200 

4 X 1200 

12 X 400 

0'.'8 

0'.'7 

(y.'9 

2011-Apr-03 

(18), (19) 

GRB 

100418A 

0.063 

4 X 1200 

4 X 1200 

8x600 

T.'O 

0'.'9 

(y.'9 

2010-Apr-21 

(20), (21) 

GRB 

100424A 

0.029 

8x600 

8x600 

8x600 

T.'O 

0'.'9 

(y.'9 

2013-Mar-ll 


GRB 

100508A 

0.024 

2 X 1200 

2 X 1200 

8 x300 

I'.'O 

0'.'9 

(y:9JH 

2012-Dec-25 


GRB 

100606A 

0.023 

4x900 

4x900 

9x300 

I'.'O 

0'.'9 

(y:9 

2014-NOV-18 





4x900 

4x900 

9x300 

I'.'O 

0'.'9 

(Y.'9 

2014-NOV-28 


GRB 

100615A 

0.039 

4 X 1200 

4 X 1200 

16 X 300 

T.'O 

0'.'9 

(y:9JH 

2013-Mar-05 


GRB 

100621A 

0.027 

2 X 1200 

2 X 1200 

8x300 

T.'O 

0'.'9 

a:9JH 

2012-Oct-16 

(9), (22) 

GRB 

100724A 

0.037 

4 X 1200 

4 X 1200 

12 X 400 

0'.'8 

0'.'7 

Oi'O 

2011-Apr-03 


GRB 

100728A 

0.151 

2 X 1200 

2 X 1200 

8 x300 

T.'O 

0'.'9 

(y.'9 

2012-NOV-15 


GRB 

100814A 

0.017 

4 X 1200 

4 X 1200 

8 x600 

T.'O 

0'.'9 

(Y.'9 

2010-Aug-18 


GRB 

100816A 

0.074 

4 X 1200 

4 X 1200 

8 x600 

T.'O 

0'.'9 

(Y.'9 

2010-Aug-17 


GRB 

110808A 

0.009 

4x600 

4x600 

4x600 

T.'O 

0'.'9 

(Y.'9 

2011-Aug-08 


GRB 

110818A 

0.030 

4x600 

4x600 

4x600 

T.'O 

0'.'9 

0r.'9 

2013-May-08 





4x600 

4x600 

4x600 

T.'O 

0'.'9 

(Y.'9 

2013-May-10 





8 x600 

8 x600 

8 x600 

T.'O 

0'.'9 

(y.'9 

2013-Jul-04 


GRB 

110918A 

0.017 

2 X 1200 

2 X 1200 

4x600 

T.'6 

T.'5 

(y:9JH 

2012-Dec-17 

(23) 




2 X 1200 

2 X 1200 

4x600 

T.'6 

T.'5 

(y!9JH 

2013-Jan-07 





2 X 1200 

2 X 1200 

4x600 

T.'6 

T.'5 

(y:9JH 

2013-Jan-16 


GRB 

111123A 

0.047 

4x600 

4x600 

4x600 

T.'O 

0'.'9 

Oi'O 

2013-Apr-07 


GRB 

111129A 

0.036 

4 X 1200 

4 X 1200 

16 X 300 

T.'O 

0'.'9 

Oi'O 

2011-Apr-29 


GRB 

111209A 

0.015 

4 X 2400 

4 X 2400 

16 X 600 

T.'O 

0'.'9 

(y.'9 

2011-Dec-29 

(24) 

GRB 

111211A 

0.041 

4x600 

4x600 

4x600 

T.'O 

0'.'9 

(Y.'9 

2011-Dec-13 


GRB 

111228A 

0.028 

4x600 

4x600 

4x600 

T.'O 

0'.'9 

(Y.'9 

2011-Dec-29 


GRB 

120118B 

0.124 

6x600 

6x600 

6x600 

T.'O 

T.'O 

(y:9JH 

2013-Feb-13 


GRB 

120119A 

0.096 

6x600 

6x600 

6x600 

T.'O 

T.'O 

a:6JH 

2013-Feb-26 


GRB 

120211A 

0.026 

4 X 1200 

4 X 1200 

16 X 300 

T.'O 

0'.'9 

(Y.'9 

2013-Feb-17 


GRB 

120224A 

0.032 

4x600 

4x600 

4x600 

T.'O 

0'.'9 

(y.'9 

2012-Feb-25 


GRB 

120422A 

0.030 

4 X 1200 

4 X 1200 

16 X 300 

T.'O 

0'.'9 

(Y.'9 

2012-Apr-22 

(25) 

GRB 

120624B 

0.048 

1 X 1200 

1 X 1200 

4x300 

T.'O 

0'.'9 

(Y.'9 

2012-Jul-12 

(26) 




3 X 1200 

3 X 1200 

12 X 300 

T.'O 

0'.'9 

(Y.'9 

2012-Jul-13 


GRB 

120714B 

0.008 

4 X 1200 

4 X 1200 

16 X 300 

T.'O 

0'.'9 

(Y.'9 

2012-Jul-15 


GRB 

120722A 

0.042 

2 X 2400 

2 X 2400 

8x600 

T.'O 

0'.'9 

(Y.'9 

2012-Jul-23 


GRB 

120805A 

0.027 

4x900 

4x900 

4x900 

T.'O 

T.'O 

(y:9JH 

2012-Sep-15 


GRB 

120815A 

0.099 

4x600 

4x600 

4x600 

T.'O 

0'.'9 

(Y.'9 

2013-Apr-14 

(27) 




4x600 

4x600 

4x600 

T.'O 

0'.'9 

(Y.'9 

2013-Apr-16 





4x600 

4x600 

4x600 

T.'O 

0'.'9 

(Y.'9 

2013-Apr-21 


GRB 

121024A 

0.088 

4x600 

4x600 

4x600 

T.'O 

0'.'9 

0r.'9 

2012-Oct-24 

(28) 

GRB 

121027A 

0.017 

14 X 600 

14 X 600 

14 X 600 

T.'O 

0'.'9 

(Y.'9 

2012-NOV-30 

(23) 

GRB 

121201A 

0.008 

4 X 1200 

4 X 1200 

16 X 300 

T.'O 

0'.'9 

(y:9JH 

2012-Dec-02 


GRB 

121209A 

0.039 

2x600 

2x600 

2x600 

T.'O 

0'.'9 

(Y.'9 

2012-NOV-13 


GRB 

130131B 

0.029 

8 x900 

8 x900 

24 X 300 

T.'O 

0'.'9 

(y:9JH 

2013-Mar-09 


GRB 

130427A 

0.017 

2x600 

2x600 

2x600 

T.'O 

0'.'9 

(y:9JH 

2013-Apr-28 

(29) 

GRB 

130701A 

0.073 

2x600 

2x600 

2x600 

T.'O 

0'.'9 

a.'9JH 

2013-Aug-01 


GRB 

130925A 

0.018 

4 X 1470 

4 X 1500 

20 X 300 

T.'O 

0'.'9 

a:9JH 

2013-Oct-25 

(30), (31) 

GRB 

131103A 

0.009 

4x600 

4x600 

4x600 

T.'O 

0'.'9 

(y:9JH 

2013-NOV-05 


GRB 

131105A 

0.030 

8 x600 

8 x600 

8 x600 

T.'O 

0'.'9 

(Y.'9 

2013-NOV-05 


GRB 

131231A 

0.022 

4x600 

4x600 

4x600 

T.'O 

0'.'9 

a:9JH 

2014-Feb-01 


GRB 

140114A 

0.014 

6x900 

6x900 

18 X 300 

T.'O 

0'.'9 

(y:9JH 

2014-Mar-28 


GRB 

140213A 

0.131 

2x600 

2x600 

2x600 

T.'O 

0'.'9 

a.'9JH 

2014-Feb-14 


GRB 

140301A 

0.026 

12 X 600 

12 X 600 

12 X 600 

T.'O 

0'.'9 

a:9JH 

2014-Mar-02 


GRB 

140430A 

0.117 

2x600 

2x600 

2x600 

T.'O 

0'.'9 

(y.'9 

2014-Apr-30 


GRB 

140506A 

0.082 

8 x600 

8 x600 

8 x600 

T.'O 

0'.'9 

(Y.'9 

2014-May-07 

(32) 




8 x600 

8 x600 

8 x600 

T.'O 

0'.'9 

(y.'9 

2014-Jun-08 
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Table 2: Fluxes of Balmer lines for GRB hosts 


GRB host 

Redshift 

Hd 

Hr 

Hy6 

Ha 

GRB050416A 

0.6542 

0.3+ 0.3+ 0.1 

1.4+ 0.3+ 0.2 

2.9 + 0.4 + 0.3 

16.0 + 2.0+1.6 

GRB050525A 

0.6063 

0.1 +0.1 +0.1 

0.1 +0.1 +0.1 

0.3+ 0.2+ 0.1 

0.7 + 0.2 + 0.3 

GRB050714B 

2.4383 


0.6 + 0.3 + 0.2 

0.5 + 0.2 + 0.2 

3.4+ 0.6+ 0.8 

GRB050819A 

2.5042 


0.6+ 0.3+ 0.1 

1.0+ 0.3+ 0.2 


GRB050824 

0.8277 

0.7+ 0.1 +0.1 

1.1 +0.1 +0.2 

2.6 + 0.4 + 0.5 

6.8+ 0.6+ 1.3 

GRB050915A 

2.5275 


0.5+ 0.2+ 0.1 

1.8+ 0.2+ 0.4 


GRB051001 

2.4295 


1.3+ 0.6+ 0.3 

1.9+ 0.3+ 0.4 

13+4 + 2 

GRB051016B 

0.9358 

2.4+ 0.3+ 0.2 

6.3 + 0.4 + 0.5 

10.9+ 1.2+ 1.1 

40 + 6 + 4 

GRB051022A 

0.8061 

2.8+ 0.5+ 0.3 

8.1 +0.8+ 0.7 

22.1 + 1.3+ 1.7 

110 + 4+ 15 

GRB051117B 

0.4805 

-1.4+ 0.8+ 0.1 

1.1 + l.O + O.l 

3.0+ 1.0+ 0.2 

21.1+4.1 + 1.3 

GRB060204B 

2.3393 



4.1 + 1.1 +0.6 

16.9+ 1.4+ 2.0 

GRB060306 

1.5597 



1.4+ 2.0+ 0.4 

8.9+ 3.4+ 1.4 

GRB060604 

2.1355 



0.4+ 0.1 +0.1 

1.8+ 0.3+ 0.5 

GRB060707 

3.4246 



1.1 + 1.8+ 0.8 


GRB060719 

1.5318 


0.1 +0.5+ 0.1 

0.8+ 0.5+ 0.2 

3.8+ 0.4+ 0.7 

GRB060729 

0.5429 


0.0+ 0.3+ 0.1 

0.6+ 0.2+ 0.1 

3.7+ 1.7+ 0.5 

GRB060805A 

2.3633 


0.5 + 0.4 + 0.2 

1.4+ 0.3+ 0.4 

3.7+ 0.5+ 0.9 

GRB060814 

1.9223 



8.3+ 3.0+ 1.5 

28+4 + 4 

GRB060912A 

0.9362 

0.5+ 0.3+ 0.1 

2.5 + 0.3 + 0.2 

5.9+ 0.5+ 0.5 

15.8 + 3.8 + 1.6 

GRB060923B 

1.5094 



-1.6+ 0.6+ 0.1 

2.6+ 0.8+ 0.3 

GRB060926 

3.2090 

0.1 +0.1 +0.1 


0.3 + 0.2 + 0.2 


GRB061021 

0.3453 

0.2+ 0.2+ 0.1 

0.2+ 0.2+ 0.1 

0.5 +0.1 +0.1 

1.9+ 0.1 +0.2 

GRB061110A 

0.7578 

0.0+ 0.1 +0.1 


0.3 +0.1 +0.1 


GRB061202 

2.2543 



1.2+ 0.4+ 0.3 

6.1 +0.8+ 1.3 

GRB070103 

2.6208 

1.1 + 1.9+ 0.6 

3.9+ 1.6+ 1.2 

4.3+ 0.8 + 1.1 


GRB070110 

2.3523 


0.6 + 0.4 + 0.2 


3.5+ 0.6+ 0.8 

GRB070129 

2.3384 


0.7 + 0.6 + 0.2 

1.5 +0.9+ 0.4 

6.5 + 1.1 + 1.1 

GRB070224 

1.9922 


0.1 +0.1 +0.1 

0.2+ 0.1 +0.2 


GRB070306 

1.4965 

1.7+ 1.3+ 0.1 

7.7 + 3.6 + 0.8 

11.6+ 1.2+ 0.7 

53.5+1.4 + 3.7 

GRB070318 

0.8401 


0.3+ 0.1 +0.1 

1.2+ 0.2+ 0.2 

3.3+ 0.4+ 0.5 

GRB070328 

2.0627 

0.3+ 0.2+ 0.1 

0.1 +0.3+ 0.1 

1.1 +0.2+ 0.2 


GRB070419B 

1.9586 



0.9+ 0.3+ 0.3 

4.4+ 0.6+ 1.1 

GRB070521 

2.0865 



-0.2 + 0.9 + 0.2 

10 + 5 + 3 

GRB070802 

2.4538 


0.4+ 0.2+ 0.1 

1.4+ 0.2+ 0.3 

5.1 +0.7+ 1.0 

GRB071021 

2.4515 


0.7 + 0.5 + 0.2 

2.8+ 0.3+ 0.5 

8.9+ 1.8+ 1.8 

GRB080207 

2.0856 


1.0+ 0.6+ 0.3 

1.3+ 0.5+ 0.3 

10.9+ 1.9+ 2.2 

GRB080413B 

1.1012 


0.2+ 0.1 +0.1 


2.6+ 1.0+ 0.8 

GRB080602 

1.8204 

1.6+ 0.9+ 0.3 

2.1 +0.7+ 0.3 


30+8 + 4 

GRB080605 

1.6408 



7.7+ 0.8+ 1.3 

29.1 + 1.4+ 4.3 

GRB080804 

2.2059 



0.9 + 0.4 + 0.3 

3.7+ 0.8+ 1.2 

GRB080805 

1.5052 



1.7+ 0.6+ 0.5 

11.0+ 1.6+ 2.7 

GRB081109 

0.9785 


1.8+ 0.4+ 0.3 

5.5+ 0.5+ 0.8 

20.9 + 0.9 + 2.7 

GRB081210 

2.0631 

0.3+ 0.3+ 0.1 

1.1 +0.4+ 0.3 

2.1 +0.3+ 0.5 


GRB081221 

2.2590 



2.4+ 1.3 + 1.7 

9.5 + 1.9+ 5.3 

GRB090113 

1.7494 


3.6+ 1.2+ 0.7 

3.3 +0.9+ 0.6 

15 + 2 + 3 

GRB090201 

2.1000 


3.8+ 1.8+ 1.1 

6.0+ 1.4+ 1.4 

23 + 2 + 4 

GRB090323 

3.5832 


-0.3+ 0.4+ 0.1 

0.0+ 0.5+ 0.1 


GRB090407 

1.4478 

-0.2+ 0.3+ 0.1 

-0.1 +0.3+ 0.1 

0.7+ 0.2+ 0.1 

4.4 + 0.4 + 0.6 

GRB090926B 

1.2427 

0.5+ 0.2+ 0.1 

0.7 + 0.4 + 0.2 

2.7 + 0.9 + 0.5 

14.2 + 2.5 + 2.4 

GRB091018 

0.9710 


0.0+ 0.2+ 0.1 

1.5+ 0.8+ 0.1 

4.2 + 0.6 + 0.3 

GRB091127 

0.4904 

0.7+ 0.2+ 0.1 

0.7+ 0.2+ 0.1 

1.6+ 0.2+ 0.1 

5.5+ 0.3+ 0.5 

GRB100316D 

0.0592 

14.0 + 0.7 

27.8+ 1.4 

62 + 3 

212+ 11 

GRB100418A 

0.6235 

1.6+ 1.3+ 0.3 

4.5 + 0.5 + 0.5 

10.0 + 0.7+ 1.0 

35 + 3+4 

GRB 100424A 

2.4656 


0.8 + 0.4 + 0.3 

2.0 + 0.2 + 0.5 

6.5 + 1.6+ 1.9 

GRB 100508A 

0.5201 


1.0+ 0.9+ 0.1 

7.2 + 0.6 + 0.5 

24.9+ 1.8 + 1.6 

GRB 100606A 

1.5545 



1.7+ 1.0+ 0.7 

4.7+ 0.5+ 1.4 

GRB100615A 

1.3978 


0.2+ 0.4+ 0.1 

1.0+ 0.4+ 0.2 

4.8+ 0.5+ 0.7 

GRB 100621A 

0.5426 


18.8 + 0.8+ 1.1 

43.8+ 1.0+ 2.4 

128 + 5 + 7 

GRB 100724A 

1.2890 


0.2+ 0.3+ 0.1 

1.1 +0.7+ 0.5 

3.8+ 0.5+ 1.1 

GRB 100728A 

1.5670 


1.2+ 2.5+ 1.4 

2.5+ 2.5 + 1.9 

10.6+ 1.7+ 4.7 

GRB 100814A 

1.4392 



1.3+ 0.3+ 0.2 

4.0 + 0.3 + 0.4 
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Table 2: Balmer lines fluxes (continued) 

GRB host 

Redshift 

Hb 

Hr 

n/3 

Ha 

GRB100816A 

0.8048 



1.7+ 0.4+ 0.2 

19 + 2 + 3 

GRB 110808A 

1.3490 



2.0 + 0.6 + 0.4 

7.6 + 0.5+ 1.1 

GRB110818A 

3.3609 



1.6 + 0.4 + 0.3 


GRB110918A 

0.9843 


-0.1 + 1.6 + 0.1 

10.3 + 2.4 + 0.6 

42+11+4 

GRB111123A 

3.1513 





GRB111129A 

1.0796 


1.1 +0.5+ 0.1 

3.3 + 1.9+ 0.4 


GRB111209A 

0.6770 

0.2+ 0.1 +0.1 

0.3+ 0.1 +0.1 

0.7+ 0.1 +0.1 

2.4 + 0.6 + 0.5 

GRB111211A 

0.4786 


0.5 + 0.7+ 0.1 

-0.5+ 0.4 + 0.1 

2.3 + 0.4 + 0.2 

GRB111228A 

0.7164 


0.1 +0.2 + 0.1 

0.5+0.3 + 0.1 


GRB 120118B 

2.9428 


2.0 + 0.3 + 0.3 

2.4 + 0.8 + 0.3 


GRB120119A 

1.7291 

0.2+ 1.6 + 0.3 

1.5 + 0.7 + 0.4 

4.9+ 1.0+ 1.1 

19.2+ 1.1 +3.7 

GRB 120422A 

0.2826 

3.2 + 0.3 + 0.1 

6.4 + 0.5 + 0.3 

14.8 + 0.7 + 0.6 

57.5+ 1.1 +2.1 

GRB120624B 

2.1974 

-0.2 + 0.6 + 0.1 


2.9+ 1.4+ 1.0 

10.3 + 0.7 + 3.1 

GRB 120714B 

0.3985 

0.7 + 0.5 + 0.1 

0.5 + 0.3 + 0.1 

2.5 + 0.2 + 0.2 

7.5+ 0.3+ 0.5 

GRB 120722A 

0.9590 

1.4+ 0.2 + 0.2 

2.6 + 0.2 + 0.3 

8.7 + 1.2+ 1.0 

32.0+1.5 + 3.4 

GRB120815A 

2.3587 



0.3 +0.1 +0.1 

0.8 + 0.2 + 0.2 

GRB 121024A 

2.3012 



7.3 + 0.6 + 0.6 

17.7+ 1.2+ 1.5 

GRB121027A 

1.7732 





GRB121201A 

3.3830 

0.6+ 0.4 + 0.1 




GRB130131B 

2.5393 



0.5 + 0.4 + 0.2 


GRB 130427A 

0.3401 

-0.7 + 1.0+ 0.1 

0.3 + 1.1 +0.1 

4.8 + 0.9 + 0.4 

14.6+ 1.3 + 1.1 

GRB130701A 

1.1548 



1.3 +2.4+ 0.2 


GRB 130925A 

0.3483 

1.2+ 1.0+ 0.1 

5.9+ 1.0+ 0.3 

11.7 + 0.7 + 0.6 

54 + 2 + 3 

GRB131103A 

0.5960 

3.1 + 1.1 +0.6 

6.0+ 1.3+ 1.0 

20.2+ 1.1 +2.8 

51 + 5 + 8 

GRB131105A 

1.6854 

1.0 + 0.6 + 0.3 

0.7 + 0.6 + 0.3 

1.7 + 0.5 + 0.5 

10.1+0.7 + 2.3 

GRB131231A 

0.6427 

1.1 +0.3 + 0.1 

2.4 + 0.3 + 0.2 

3.9 + 0.5 + 0.4 

14.1 + 1.1 + 1.4 

GRB140213A 

1.2079 



0.6 + 0.7 + 0.1 

1.2 + 0.5 + 0.1 

GRB140301A 

1.4155 

2.1 +0.6 + 0.2 

1.4 + 0.6 + 0.2 

5.0 + 0.6 + 0.5 

31.7+1.5+2.7 

GRB140430A 

1.6019 




5.8+ 0.8+ 0.6 

GRB140506A 

0.8893 




1.2+ 0.4+ 0.2 


Notes. Measurements are in units of 10“'’ ergcm“’ s“', and are corrected for the Galactic foreground reddening. A 
correction for slit-loss based on broad-band photometry has been applied to the measurements as described in the 
text. These values thus represent host-integrated measurements except for GRB 100316D (see below). Line fluxes are 
not corrected for host-intrinsic extinction. The first error represents the statistical error due to photon statistics and 
line-flux measurement. The second error is the systematic error in the absolute flux calibration due to slit-loss and 
scaling to photometry. Redshifts are given in a heliocentric reference frame. No data means that either the wavelength 
range of the respective line is not covered, all data in that wavelength range have been excluded from the automated 
fitting procedure, or no meaningful constraints could be obtained for the given line. 

" GRB 100316D has the lowest redshift in the sample (z - 0.0592). The fraction of the host that is covered by the slit 
is too small to derive representative host-integrated measurements. These values are thus as derived from the observed 
spectrum and not scaled by photometry. No systematic error on the line fluxes is given in this case. 
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Table 3: Fluxes of forbidden lines (continued) 
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Table 4: Physical properties of GRB hosts 


GRB host'* 

Redshift 

Eb-v 

SF-tracer 

SFR* 

cr 

Z 



(mag) 


(Moyr-‘) 

(kms“‘) 

12 + log(0/H) 

GRB050416A 

0.6542 

““q^q+UTT" 

Ha 

4.5!|;° 

47 + 4 

8.46!°;|| 

GRB050525A 

0.6063 

^■^^- 0.10 

Ha 

-0.05 

26 + 5 


GRB050714B 

2.4383 

0 21 +“-^^ 
'^■^^- 0.21 

Ha 

19 q+14.0 

35 + 16 


GRB050819A 

2.5042 

A ^4+0.84 
-0.34 

H/1 

99+426 

-15 

50 + 6 


GRB050824 

0.8277 


Ha 

120 +°™ 
^■^'+-0.26 

48 + 5 

Oil +0.18 
- 0.20 

GRB050915A 

2.5275 


W 

196+1563 

^^^-174 

87 + 10 


GRB051001 

2.4295 

0.58!“« 

Ha 

110+124 
iiU 55 

67 + 8 


GRB051016B 

0.9358 

0 05 

Ha 

1 0 9 + 2.6 
- 2.0 

58 + 4 

0 97+0.15 
- 0.20 

GRB051022A 

0.8061 

0 

^•'^^-0.04 

Ha 

60 +™ 

88 + 5 

0 49+0.09 

o-^^_0.09 

GRB051117B 

0.4805 

0 72+°-^’ 

'^■'^-0.24 

Ha 

4 7 + 4.9 
^•'-2.2 

85 + 11 

9.00!°;™ 

GRB060204B 

2.3393 

0 34+0-29 
-0.23 

Ha 

78!«5 

85 + 8 


GRB060306 

1.5597 

0-44!°:® 

Ha 

17.6!*?;° 

60+14 


GRB060604 

2.1355 

0.38!“2 

Ha 

7 9+9.4 
'•^-3.6 

68 + 13 

8.10!°;/ 

GRB060707 

3.4246 


[On] 

IQ Q+48.0 
^^•^-14.3 



GRB060719 

1.5318 

^ -0.36 

Ha 

7 1 +18.9 
'•^-3.9 

42 + 5 

8.61!°;/ 

GRB060729 

0.5429 

0 71+0-43 
^■^^-0.47 

Ha 

0 96+2.21 
^■^^-0.69 

66 + 16 


GRB060805A 

2.3633 

^■^^- 0.00 

Ha 

9 o+^° 

^■'+- 2.5 

71+7 


GRB060814 

1.9223 

0 17+0-39 
U.i /_o.i7 

Ha 

54+89 

-19 

132+ 11 

8.61!°;/ 

GRB060912A 

0.9362 

-0.09 

Ha 


62 + 5 

GRB060923B 

1.5094 


Ha 

*^•^-1.5 

46 + 15 


GRB060926 

3.2090 


[Om] 

26+"^"^ 

4.0_17 

60 + 8 


GRB061021 

0.3453 

0 -lltn 

Ha 

0 05 

^■^-’-001 

16.1+4.8 


GRB061110A 

0.7578 


H/3 

0 23+0.38 
-0.15 

31+4 


GRB061202 

2.2543 

0.58!“4 

Ha 

49+60 

^*^-22 

64 + 7 


GRB070103 

2.6208 

- 0.00 

H/3 

49+162 

^■^-17 

124 + 30 


GRB070110 

2.3523 

^■^^- 0.00 

Ha 

8.9!™,° 

24 + 4 


GRB070129 

2.3384 

0 17+0-35 

U.i /_0.17 

Ha 

20 !f 

76+11 


GRB070224 

1.9922 


[Om] 

9 9 + 6.5 
-2.3 

37 + 9 


GRB070306 

1.4965 

A 40 + 0.08 

^■^*^-0.07 

Ha 

101 +^"^ 

121+55 

0 54+0.09 
^•'^^-0.09 

GRB070318 

0.8401 

^■^'^- 0.14 

Ha 

A 70 + 0.44 

^■'^-0.24 

53 + 5 


GRB070328 

2.0627 

0 16+°™ 
'+■^+’-0.16 

H/3 

8.4+™°-" 

93 + 14 


GRB070419B 

1.9586 

^■-’^-0.30 

Ha 

21 +™ 

86 + 10 


GRB070521 

2.0865 


Ha 

26!?^ 

249 + 108 


GRB070802 

2.4538 

^■■^^- 0.12 

Ha 

24!™ 

57 + 5 


GRB071021 

2.4515 

0 19+046 
^■^^-0.17 

Ha 

09+20 

-12 

100 + 17 


GRB080207 

2.0856 

^■^^-0.25 

Ha 

77+86 
' '-38 

136+18 

0 74+0.15 

GRB080413B 

1.1012 

0.43!°i° 

Ha 

21 +^+ 

39 + 5 

0 29+0.32 
o-^^_0.30 

GRB080602 

1.8204 

0-58!°:™ 

Ha 

125 + 145 

91 + 13 

0 54+0.09 
^•''^-0.09 

GRB080605 

1.6408 

0-26tlo 

Ha 

47!™ 

80 + 6 

GRB080804 

2.2059 

0-38!SJ^ 

Ha 

15.2!4i/ 

50 + 9 

0 49+0.13 

GRB080805 

1.5052 

0.78!°;™ 

Ha 

45!^^ 

54+12 

GRB081109 

0.9785 

0 36+^-^^ 

Ha 

11 . 8 !^;^ 

108 + 6 

0 7 C+ 0.09 
^•'''-0.09 

GRB081210 

2.0631 

^■^■^-0.13 

H^ 


118+ 12 


GRB081221 

2.2590 

0 31+0-55 

Ha 

oc+106 
■2 *^-22 

93 + 12 


GRB090113 

1.7494 

A 02+0.22 
^•^^- 0.01 

Ha 

17 A+10.1 
^ '-^-4.8 

70 + 9 


GRB090201 

2.1000 

0 11 +°™ 

Ha 

48+3° 

^°-14 

171 + 12 


GRB090323 

3.5832 


[On] 

24+53 
^ -17 

60+13 


GRB090407 

1.4478 

-0.26 

Ha 

13.8!™/ 

109 + 8 

0 qc+0.13 
°*°^-0.13 

GRB090926B 

1.2427 

0 63+^-^^ 
^■^■^-0.18 

Ha 

^^-11 

65 + 4 

0 34+0.15 
^•■^^-0.17 

GRB091018 

0.9710 

^■^^-0.06 

Ha 

1 9Q+3.46 
^■^^-0.32 

57+10 

Q 7Q+0.18 

^•'°-0.19 
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Table 4: Physical properties of GRB hosts (continued) 

GRB host 

Redshift 

Eb-v 

SF-tracer 

SFR“ 

cr 

Z 



(mag) 


(Moyr-‘) 

(kms-‘) 

12 + log(0/H) 

GRB091127 

0.4904 

^■^^-0.08 

Ha 

A 92+0.10 

-0.07 

30 + 5 

8.07!i;;- 

GRB 100418A 

0.6235 

0 17+006 
-0.07 

Ha 

4 9+1.0 
^•^-0.8 

56 + 4 

0 C 9 +O.IO 
- 0.10 

GRB100424A 

2.4656 

0 13+018 

H^ 

2 i:f 

87 + 5 

7 07+0-25 
^•^*^-0.18 

GRB100508A 

0.5201 

A 9Q+0.09 
^■^^-0.09 

Ha 

2 - 6 !°:’ 

80+11 

8 . 68 !?;!° 

GRB 100606A 

1.5545 

0-05!°:“ 

Ha 

4 q+12.9 

107 + 36 

8.71!°;?? 

GRB 100615A 

1.3978 

0.48!“« 

Ha 

8 . 6 !^” 

45 + 5 

8.40!?;!? 

GRB100621A 

0.5426 

0 05 

Ha 

00 

1 + 
o c 
bo be 

82 + 4 

0 C 9 -FO.IO 

GRB100724A 

1.2890 

-0.24 

Ha 

3 2+^-* 

58 + 7 


GRB100728A 

1.5670 

A 99+0.69 
^■^*^-0.23 

Ha 

14 c+60.6 

57 + 8 


GRB 100814A 

1.4392 

O-O^loM 

Ha 

3 2+^-^ 

-0.7 

31+5 


GRB100816A 

0.8048 

1 99+0.24 
^■■^^-0.22 

Ha 

58^^ 

111+30 

Q '7C+0.16 
^•'^-0.18 

GRB110808A 

1.3490 

-0.25 

Ha 

8-3!1'6' 

40 + 4 

7 09+0.31 
^•^*^-0.23 

GRB110818A 

3.3609 

^■'^'^-0.31 

H/3 

44!® 

89 + 8 

Q 9C+O.I7 
°*^^-0.25 

GRB110918A 

0.9843 

Ha 

99+28 

126+ 18 

8.93!?;!! 

GRB111123A 

3.1513 


[On] 

77+163 
' -52 

135 + 21 

8 - 01 !?:?? 

GRB111129A 

1.0796 


[On] 

^ 1+10.8 
'^•■^-3.4 

117 + 35 

GRB 111209A 

0.6770 


Ha 

^■■^'^-0.13 

35 + 5 

7 ac+0.30 

17 

GRB111211A 

0.4786 

^■^^- 0.00 

Ha 

-0.03 

38 + 8 


GRB111228A 

0.7164 


H/3 

0.32!°!^ 

19.7 + 5.5 


GRB120118B 

2.9428 

^■^^- 0.00 

up 

28!?! 

193 + 8 

7.89!?;?? 

GRB 120119A 

1.7291 

^■■^'^-0.14 

Ha 

43+24 

^■7-14 

104 + 17 

8.60!?;!? 

GRB 120422A 

0.2826 

A 27+0.03 
-0.03 

Ha 

1 9Q+0.13 
^■■^^-0.12 

25 + 4 

Q 90+0.09 
^•*^^-0.09 

GRB120624B 

2.1974 

u.zi 0.21 

Ha 

30!?? 

77 + 6 

0 49+0.20 
^•^*^-0.27 

GRB 120714B 

0.3985 

A ^a+0.08 
-0.08 

Ha 

0.27!°;® 

34 + 4 

8.39!?;!! 

GRB120722A 

0.9590 

^•^^-0.05 

Ha 

22 !^ 

56 + 4 

8.48!?;!? 

GRB 120815A 

2.3587 

^■^^-0.06 

Ha 

7 it+2.7 
^•■2-1.0 

28 + 5 

0 4 J+ 0.11 

^•^^- 0.12 

GRB121024A 

2.3012 

^■^^- 0.00 

Ha 

37+4 
■7/-4 

88 + 4 

GRB 121027A 

1.7732 


[Om] 

24+41 

z.His 

119 + 75 


GRB121201A 

3.3830 


[On] 

30!® 

86 + 17 


GRB130131B 

2.5393 


[Om] 

8 - 0 !??," 

73+29 


GRB130427A 

0.3401 

^■^^-0.06 

Ha 

0 - 34 !°:?° 

40 + 5 

Q C7+O.I2 
-0.13 

GRB130701A 

1.1548 


[On] 

0 - 78 !?:® 

82 + 42 


GRB 130925A 

0.3483 

0 41+006 

Ha 

o q+0.5 

49 + 5 

0 79+0.08 
^•^■^-0.08 

GRB131103A 

0.5960 

^■^^-0.06 

Ha 

4 4+l‘^ 
^■^-0.9 

87 + 7 

8.48!?;!? 

GRB131105A 

1.6854 

^•'^'^-0.18 

Ha 

■^^-13 

52+11 

8.61!?;?? 

GRB131231A 

0.6427 

^■^^-0.02 

Ha 

1 90+U.28 
■■70-0 20 

33+4 

8 

^'^'^- 0.12 

GRB 140213A 

1.2079 

^■^^-0.06 

Ha 

0.12:111 

34+14 


GRB 140301A 

1.4155 

''^-O.lO 

Ha 

106 !?? 

117 + 6 

8.89!?;?? 

GRB140430A 

1.6019 


Ha 

0 C+7.I 
°-^-3.8 

40 + 7 

8 - 67 !?:!? 

GRB140506A 

0.8893 


Ha 

^■'7^-019 

61+9 



Notes. “ The physical parameters presented here are integrated and thus averaged over the entire galaxy. 
We do not perform a resolved analysis (neither spatially, nor in velocity space). In some cases (e.g., 
GRB 120422A, Schulze et al.|2014 i, a spatially resolved analysis leads to somewhat different results and 
interpretation of the galaxy properties. * The quoted error on SFR is logarithmic because it contains the 
error in the dust correction. The derived SFR also has a lower limit because of the physical condition that 
Ay > Omag. This lower limit is given by SFR^j^ = 4.8 x FHt,, 42 - 
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Appendix A: Notes on individual targets 

Appendix A. 1: GRB 111129A 

The X-Shooter spectrum of the GRB 111129A afterglow 
( Racusin et al.|20lT| l was taken on 2011-Nov-30, approximately 
8.3 h after the initial Swift trigger. Above a red continuum ( |Rossi| 
|et al.||20TT] l we detect a single emission line (significance of 
5.4cr) at a wavelength of 7756 A. Interpreting this as [On], the 
redshift of GRB 111129A is z = 1.0796. If it were any of the 
other strong lines, we would expect to detect at least one other 
emission line in the spectrum, which we do not. At z = 1.0796, 
however, there are indications for [O iii] and H/3 at the 2cr level 
in the NIR arm (Ha is in the 7//-bandgap and thus not seen). 


Appendix A.2: GRB 120211A 

In deep Keck imaging taken on 2013-Feb-10, we detect an 
/ab = 25.0 + 0.2maghostgalaxy candidate within the XRTerror 
circle dGoad et al.|20T2l l of GRB 120211A dSonbas et al.|20T2l l. 
The spectral continuum of the putative host galaxy is detected 
also at high significance in the X-Shooter spectrum from 2013- 
Mar-20 in a wavelength range between 9000 A and 4100 A. Be¬ 
low 4100 A the flux sharply drops to zero. This apparent break 
is consistent with the onset of Lya absorption at z ~ 2.4. Con¬ 
sistent with this redshift, there are minor indications at the 2cr- 
significance level for [Oiii]/I5007 emission at z = 2.346 located 
in a region of high sky emission. 


Appendix A.3: GRB 120224A 


In the X-Shooter spectrum ( Wiersema et al.|2012b i of an after¬ 
glow candidate of GRB 120224A ( Saxton et al.|2012] l taken on 
2012-Feb-25, we detect a red continuum in the VIS and NIR 
arms, without obvious emission lines. Later GROND imaging 
reveals this to be a very red (r - Kab =3.7 mag), and constant 
source, a plausible GRB host candidate. The spectral continuum 


drops below the noise floor blueward of 7000 A. We rule out that 
the break is due to the Lya at z = 4.8 because the photometric 
colors, in particular the ^-band detection (g' = 25.2 + 0.2) disfa¬ 
vors such a high redshift. A photometric redshift analysis gives 
90% confidence interval of 0.9 < Zphot < 1-3. Indeed, there is a 
small redshift window at z ~ 1.1 in which all strong emission 
lines ([On], [Oiii] and Ha) would be hidden in telluric bands. 
This would provide an explanation for the non-detection of emis¬ 
sion lines in our spectrum. The drop to zero flux in the spectrum 
would in this case be caused by intrinsic redness of the galaxy 
and the Balmer break. Significantly higher or lower redshifts are 
unlikely: The stellar mass of the galaxy derived from the pho¬ 
tometry exceeds 10**Mo already at z > 1.5, making the non¬ 
detection of emission lines hard to explain. Similarly, at z < 0.6, 
we would have expected to detect at least Ha in the spectrum 
even for a dusty galaxy. 


Appendix A.4: GRB 120805A 


The galaxy coincident with the afterglow (Guidorzi & Mundell 

2012 

[Gorosabel et al.||20121l of GRB 120805A (Troja et al. 

2012 

is identified through NOT (|Malesani et al. |2012|| and 


late GROND and VLT/HAWK-I observations at a brightness of 
r = 24.1 ± 0.1 mag. It was observed with X-Shooter spectro¬ 
scopically on 2012-Aug-15. The spectral continuum is detected 
at high significance in the UVB/VIS arm between 5000 A and 


9500 A. No signal is recorded blueward of 5000 A or in the 
NIR arms. The shape of the spectral continuum and the GROND 
photometry are consistent with a redshift of z ~ 3.1. The Lya 
break explains the drop in flux at around 5000 A, and the ab¬ 
sence of strong emission lines in the spectrum is caused by tel¬ 
luric absorption bands. At this redshift, the host is luminous: 

M| 7 nnA-21.9 mag corresponds to ~ 2.5 ■ M* at z ~ 3 (e.g., 

Reddy & Steidel|2009| l. 


Appendix A.5: GRB 121209A 


The Rab ^ 24.1 + 0.1 mag host galaxy of GRB 12 1209A 
([Maselli et al.|20T^ was ident ified already very early ( |Kruh-| 
|ler et ak 2012[ |Perley et al.||2012| l, and observed with X- 
Shooter on 2012-Dec-13. The galaxy-continuum is detected sig¬ 
nificantly with X-Shooter between 4900 A and 9500 A but not 
redward and with Keck LRIS b etween the atmospheric cutoff 
and 10300 AlPerley et al.||2012i. Using all available late pho¬ 
tometry from VLT/FORS2 and VLT/HAWK-I, the best-fit pho¬ 
tometric redshift for this galaxy i s Zphot = 2.1+ 0.3, in agreement 
with earlier limits ( [Perley et al.||2()12 i. Despite the bright con¬ 
tinuum, no emission lines are seen. This is also consistent with 
Zphot ~ 1.9 because all strong emission lines ([Out], [On], and 
Ha) would be located in telluric absorption bands at this redshift. 


Appendix A.6: GRB 140114A 


At the position of the optical af terglow ([Butler et al.|2014] |Cano| 
et al.||2014| l of GRB 140114A ( |Troja et al.||2014 l, we detect an 


Rab = 24.4 + 0.2 mag galaxy with NOT/ALFOSC that we in¬ 
terpret as the GRB’s host galaxy. Our X-Shooter spectroscopy 
from 2014-Mar-28 reveals no emission lines, but the galaxy con¬ 


tinuum between 4850 A and 9800 A. A sharp drop in flux blue¬ 
ward of this wavelength range is best described by the Lya break 
at z = 3.0. The relatively high redshift is also consistent with the 
absence of emission lines because the strong forbidden lines of 
[O It] and [O iii] are in telluric absorption bands and Ha outside 
of the wavelength response of X-Shooter. 


Appendix B: Details on metallicity measurements 
and notable individual events 

To derive oxygen abundances via the strong-line diagnostic ra¬ 
tios from Nagao et al. (20061 and Maiolino et al.| p008| l, we 
take into account the dust attenuation, its uncertainty ([N ii] ver¬ 
sus Ha or [Ne iii] versus [O ii] have a negligible dependence 
on Ay) and the systematic scatter in the respective relation. We 
use [N ii]/Ha if available to discriminate between the upper and 
lower R 23 branch. Table [^provides the calculated metallicities 
and errors (as the range between 16% and 84% of cumulative 
probability distribution). These oxygen abundances need to be 
taken with a grain of salt, in particular at higher redshifts, as in¬ 
dicated in the main text. 

In some cases, in particular at the highest redshifts (z > 2.5), 
dust corrections to the line fluxes are not available because of 
lacking Balmer lines. In those cases, however, metallicities are 
derived using line ratios that are only marginally sensitive to dust 
corrections such as [Nein]over [On]. The missing constraints 
on Eb-v are thus not a strong concern for the metallicity mea¬ 
surement. The systematic and statistical error of the [Ne iii] over 
[O It] ratio is large, and it does not provide strong additional con¬ 
straints once [N ii] -based line diagnostics are available. It is, 
however, the only accessible method at z > 3 so we use it with 
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the caveat that it is the least tested and used strong-line diag¬ 
nostic ratio. It did provide a consistent metallicity in those cases 
where we could test it with other methods. 

Figure [BTT] shows several representative examples of metal¬ 
licity determinations for six galaxies at different redshifts. We 
specifically pick GRB hosts or afterglows that were exten¬ 
sively discussed in the recent literature: GRB 070306 (|Jaun- 
sen et al.||20081 |Kruhler et al.||2011|l, GRB 080207 Punt et al. 
2011[|Svensson et al.|20121l, GRB 091127 (jVergani et al.|201 If, 
GRB 100621A Preiner et al l|20T3|), GRB 111209A dGendre 


2015 


2014 


et_^ 2013t[Ovan et al.|2014bl Greiner et al.|2015bl|Kann et al. 


GRB 070306 (z = 1.497) is a well-studied dusty GRB ( Jaunsen| 
et al.||20()8)l, with a host galaxy that is detected throughout the 


electromagnetic spectrum (Hunt et al.|2014[ Schady et al.|2014 


Appendix B.3: GRB 091127 

GRB 091127 (z = 0.490) is a typical low-redshift GRB with a 
bright afterglow and faint host ( |Vergani et al.|20lT| Filgas et al. 
|2011| l. It is an example of a galaxy having R 23 in the overlap 
region between the two metallicity solutions. Even though the 
forbidden oxygen and Balmer lines are detected at good S/N, the 
lack of detection of either [Ne iii] or [N ii] limits the number of 
useful line-diagnostics to R 23 - Given that the line ratios put the 
host of GRB 091127 at around the turn-over point of the two 
branches, the probability distribution is not double valued, but 
rather broad because 12-1- log(0/H) is not very sensitive to vari¬ 
ations in R 23 in this region. We derive 12 + log(0/H) = 8.07);g2Q 
for the host of GRB 091127. 


and GRB 130427A (e.g., |Xu et al.||2013[ |Perley et al. 
. We will discuss them briefly in the following section as 
they are illustrative for the metallicity measurement procedure 
and interesting individually. 


Appendix B.1: GRB 070306 


Appendix BA: GRB 100621A 

GRB 100621A (z = 0. 543) is a low-z GR B with a faint and 
dust-reddened afterglow preiner et al.|2013 1 . The host spectrum 
has, in addition to well-detected [On], [Out], Ha, and H^S emis¬ 
sion lines, a significant detection of [Ne iii] (but not [N ii]). The 
degeneracy of the double-peaked R 23 probability distribution in 
metallicity is clearly evident in Fig. |B.1| It is, however, broken 
through [Ne iii]. Even though the individual [Ne ni]-based metal¬ 
licity is not well constrained, in combination with R 23 we can de¬ 
rive 12-i-log(0/H) accurately to 12-Hlog(0/H) = 8.52+0.10 with 
only a very small remaining probability for the lower branch so¬ 
lution of the R 23 diagnostic. 


Appendix B.5: GRB 111209A 

GRB 111209A (z - 0.677) is one of the longest GRBs ever de¬ 
tected. The extreme temporal properties were used to postulate 
a new class of GRB s pendre et al.|2013[ [Levan et al.|2014bT l, a 
result questioned by |Virgili et al. ( 2013| l. The ultra-long duration 
is possibly indicative of a different progenitor channel (Nakauchi 
et al. 2013|l, and the associated luminous supernova is well 


|Perley et al.||2015b| . The X-Shooter spectrum shows numerous 
emission lines, which allow the Balmer decrement to be reason¬ 
ably well measured (Eb-v = 0.43];gQ2 mag). The [Nii]/Haline- 
flux ratio uniquely puts the R 23 solution onto the upper branch. 
The individual emission line diagnostics have overlapping prob¬ 
ability distributions, and a simultaneous fit to all five line ratios 
results in a metallicity of 12 + log(0/H) = 8.54 + 0.09, which is 
0.8 + 0.2 times the solar value (Eig. [B^. 


Appendix B.2: GRB 080207 

GRB 080207 (z ^ 2.086) is also a we ll-discussed dark GRB 
( Hunt et al.|201 1 Svensson et al.|20T2 i. Its host is detected by 
various facilities from the optical to radio ( [Rossi et al.||2012] 
Perley et al. 2013b| l, and has a very high stellar mass (M* ~ 
10" Mq), extremely red optical/NIR colors and is vigorously 
forming stars {SFRiia ~ 90Moyr“'). Because of the large 
Ay ~ 2 , [On] is not detected and the applied strong-line diag¬ 
nostics are thus based exclusively on [N ii]. We measure an oxy¬ 
gen abundance around solar, 12 + log(0/H) = 8.74 + 0.15. This 
is the GRB host with the highest metallicity at z > 2 in the sam¬ 
ple. Its metallicity, however, is averaged over the multiple knots 
that are detected at the spatial resolution of HST ( jSvensson et al.j 
2012|l, which in theory could host a range of metallicities. 


explained in a magnetar scenario ( jKann et al.|2015| Greiner et al.j 
|2015b| l. Based on the detection of [O ii], [O iii] and [Ne iii], in ad¬ 
dition to the Balmer lines, we measure 12+log(0/H) = 7.95];° 
or between 10% and 40% of the solar value in the adopted metal¬ 
licity scale. This puts the galaxy at the low end of metallicities 
compared to many GRBs of more normal duration at a similar 
redshift (Eig.[T7)l. Together with the low reddening, narrow line 
width, compact morphology and low luminosity pevan et al. 
2014b|l it shares many similarities with Blue Compact Dwarfs 


(e.g., [Gil de Paz et al.|[2003] and references therein) or hosts 
of super-luminous supemovae (e.g., Leloudas et al.||20T5] l. The 
galaxy hosting GRB 111209A is, however, not extremely metal- 
poor, and its spectroscopic properties do not provide a stark con¬ 
trast to long GRBs with a more common duration (see, e.g., 
GRB 091127 above), a conclusion similar to the one of|L evan| 
[e^(|20^. 

Appendix B.6: GRB 130427A 

GRB 130427A (z = 0.340) is one of the brightest GRBs of 
all time due to its low redshift and is extremely well studied 
from the radio to GeV energies (e.g., Perley et al.||2014[ Ack- 


jermann et al.||2014{ [van der Horst et al.||2014 1. It is a very rare 

example of a low-redshift GRB with an energy budget compa¬ 
rable to common high-redshift GRBs and therefore elucidates 
cosmological GRBs and their connection to SNe through a lo¬ 
cal analogue ( Xu et al.|2013 Levan et al.|2014a] Melandri et al. 
|2014j l. The host characteristics are relatively typical for the given 
redshift: With a SFRyf^ - 0.3;^° j Mq yr^'and metallicity of 


u u J jvHq. - j 

12 + log(0/H) = 8.57];^ I 2 it is well within the distribution of 
host properties for low-redshift GRBs. 
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Fig. B.l: Example of metallicity measurements using different strong line diagnostics (Nagao et al.||2006 Maiolino et al.|[2008 


Every colored line represents the normalized probability distribution of one strong-line diagnostic. The black line is the probability 
distribution of 12 -H log(0/H) simultaneously minimized against all available measured line ratios. 
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